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ABSTRACT 

To assess how external factors such as local interactions and fresh gas accretion influence the global 
ISM of galaxies, we analyze the relationship between recent enhancements of central star formation and 
total molecular-to-atomic (H2/HI) gas ratios, using a broad sample of field galaxies spanning early- 
to-late type morphologies, stellar masses of lO'''^ — 10^^''^ Mq, and diverse stages of evolution. We find 
that galaxies occupy several loci in a "fueling diagram" that plots H2/III ratio vs. mass-corrected 
blue-centeredncss, a metric tracing the degree to which galaxies have bluer centers than the average 
galaxy at their stellar mass. Spiral galaxies of all stellar masses show a positive correlation between 
II2/III ratio and mass-corrected blue-centeredness. When combined with previous results linking 
mass-corrected blue-centeredness to external perturbations, this correlation suggests a systematic link 
between local galaxy interactions and molecular gas inflow/replenishment. Intriguingly, E/SO galaxies 
show a more complex picture: some follow the same correlation, some are quenched, and a distinct 
population of blue-sequence E/SO galaxies (with masses below key scales associated with transitions in 
gas richness) defines a separate loop in the fueling diagram. This population appears to be composed of 
low-mass merger remnants currently in late- or post-starburst states, in which the burst first consumes 
the H2 while the galaxy center keeps getting bluer, then exhausts the H2, at which point the burst 
population reddens as it ages. Multiple lines of evidence suggest connected evolutionary sequences in 
the fueling diagram. In particular, tracking total gas-to-stellar mass ratios within the fueling diagram 
provides evidence of fresh gas accretion onto low-mass E/SOs emerging from their central starburst 
episodes. Drawing on a comprehensive literature search, we suggest that virtually all galaxies follow 
the same evolutionary patterns found in our broad sample. 
Keywords: Galaxies: General — Galaxies: Interactions — Galaxies: ISM 
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1. INTRODUCTION 

It has been well documented that stars form in molec- 
ular gas (e.g., Bigiel et al. 2008). Therefore, understand- 
ing what drives the conversion of hydrogen between its 
atomic and molecular forms is key to understanding how 
galaxies evolve. The molecular-to-atomic gas mass ratio 
(II2/III) varies widely between galaxies, and several stud- 
ies have aimed to determine what properties of galaxies 
- or their environments - play the largest role in the 
evolution of II2/III. 

Work from the last few decades has revealed correla- 
tions between global H2/III and such properties as lu- 
minosity (or stellar mass), total gas mass, morphology, 
and specific star formation rate (Kenney & Young 1989; 
Thronson et al. 1989; Young & Knezek 1989; Braine & 
Combes 1993; Casoli et al. 1998; Boselli et al. 2002; 
Obreschkow & Rawlings 2009). Unfortunately, these 
studies have been largely focused on only certain types of 
galaxies (e.g., massive spirals) and/or star-forming/FIR- 
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bright galaxies in the nearby universe, which are not rep- 
resentative of the galaxy population as a whole. More re- 
cently, the CO Legacy Database for the GALEX Arecibo 
SDSS Survey (COLD CASS) has measured CO emission 
for a randomly selected sample of ~300 galaxies with 
stellar masses Af, > 10^°Mq. COLD GASS finds cor- 
relations between H2/III and structural properties like 
stellar mass, global stellar mass surface density, r-band 
light concentration index, and NUV — r color, while also 
finding that there are thresholds in galaxy concentration 
index and global stellar surface mass density above which 
detections of molecular gas begin to disappear (Saintonge 
et al. 2011; Kauffmann et al. 2012). While the survey's 
selection improves upon past work by focusing on more 
"normal" galaxies, a caveat is that it only includes mas- 
sive galaxies. Moreover, multiple physical mechanisms 
may underlie the observed correlations. 

The question remains as to what balance of internal 
and external processes regulates H2/HI in galaxies. A 
number of authors have aimed to explain H2/III com- 
pletely in terms of local physics within galaxies, set by 
their structure and dynamics. Detailed studies of atomic 
and molecular gas in nearby galaxies on kpc scales (Re- 
gan et al. 2001; Kuno et al. 2007; Waher et al. 2008; 
Leroy et al. 2009) have shown that H2/III is correlated 
with local internal variables such as galactocentric ra- 
dius, stellar and gas surface mass density, radially vary- 
ing velocity dispersion and rotation velocity, or combina- 
tions of these (Wong & Blitz 2002; Blitz & Rosolowsky 
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2004; Lcroy et al. 2008; Lucero & Young 2009; Schruba 
ct al. 2011). Dwarf galaxies often show deviations from 
these relations, although such deviations are commonly 
assumed to reflect the metallicity dependence of the CO- 
to-H2 conversion factor, Xco (Wilson 1995; Arimoto 
et al. 1996; Bolatto et al. 2008; Obreschkow & Rawl- 
ings 2009; Wolfire et al. 2010; Glover & Mac Low 2011; 
Leroy et al. 2011). Correlations of H2/HI with inter- 
nal properties can be used to support models where the 
molecular fraction of gas is governed by the equilibrium 
between molecule formation on dust grains and destruc- 
tion by the FUV background (Elmegreen 1993; Blitz & 
Rosolowsky 2006; Krumholz et al. 2009). Alternative 
theories argue that molecular cloud formation is a very 
non-equilibrium process, spurred by gravitational insta- 
bilities, converging gas flows, and/or cloud-cloud colli- 
sions (Tan 2000; Mac Low & Klessen 2004; Hcitsch & 
Hartmann 2008; Pelupessy & Papadopoulos 2009), and 
that newly formed molecular clouds fail to reach equilib- 
rium with their environments before they are destroyed 
by star formation (Mac Low & Glover 2012). In this pic- 
ture, molecular cloud formation may be highly dependent 
on drivers that disrupt gas equilibrium. 

Consideration of global dynamical states suggests that 
external factors lead to increases in H2. At the very 
least, external perturbations and/or bars can transport 
gas to the central regions of galaxies where it can ex- 
ist at higher densities that promote the conversion of 
HI to H2 and then stars. Barred galaxies, which may be 
linked to interactions (Gcrin et al. 1990; Miwa & Noguchi 
1998), display central gas concentrations indicative of 
inflows (Sakamoto et al. 1999; Sheth et al. 2005). In- 
teracting systems display deviations from typical H2/III 
and star formation rate relations, with higher average 
II2/III than non-interacting galaxies (Kenney & Young 
1989; Braine & Combes 1993; Lisenfeld et al. 2011) and 
higher star formation rate density than predicted using 
the Kennicutt-Schmidt relation defined by normal spirals 
(Kennicutt 1998; Bigiel et al. 2008). Moreover, Blitz & 
Rosolowsky (2006) find higher H2/III ratios at a given 
mid-plane pressure for interacting systems compared to 
non-interacting systems. Fosi-starburst galaxies simi- 
larly show excess star formation given their gas content, 
likely due to H2 heating and depletion by the starburst 
(Leroy et al. 2006; Robertson & Kravtsov 2008; Wei et al. 
2010b). Papadopoulos & Pelupessy (2010) recreate such 
deviations in simulations, linking them to quickly evolv- 
ing systems, particularly those that are gas rich and have 
experienced recent minor mergers or fresh gas infall. It 
should also be noted that Xco niay be lower for cen- 
trally concentrated molecular gas, possibly making the 
increase in molecular gas due to an interaction appear 
even higher (Downes et al. 1993; Garcia-Burillo ct al. 
1993; Regan 2000). 

The correlation of galaxy interactions with molecular 
gas enhancement (and/or higher CO luminosity per unit 
molecular gas mass) is an established result, but an open 
question remains whether interaction-induced enhance- 
ments are occasional serendipitous events or the domi- 
nant driver of observed H2/III ratios in galaxies. Kan- 
nappan et al. (2004) argue that galaxy interactions ac- 
count for the majority of recent central star formation 
enhancements that produce blue central color excesses, 
based on the correlation of blue-centered galaxies and 



signs of minor mergers/encounters, a result that is con- 
firmed by Gonzalez- Perez et al. (2011). This correlation 
implies that blue-centered galaxies experience gas inflows 
to feed these star-forming events, a scenario addition- 
ally supported by such galaxies' transient decreased cen- 
tral metallicities (Kewley et al. 2006, 2010; Rupke et al. 
2010). Kannappan et al. (2004) find that roughly 10% of 
star-forming galaxies show blue-centered color gradients, 
and most have likely experienced a blue-centered phase 
at least once in their lifetimes; this percentage increases 
with decreasing luminosity. Additionally, when the sys- 
tematic trend towards central reddening at higher lumi- 
nosity is subtracted (yielding luminosity-corrected blue- 
centeredness), the number of galaxies identified as having 
enhanced central blueness at fixed luminosity increases, 
and the correlation with interactions is stronger. These 
results are interesting when combined with recent anal- 
yses of the rates of nfinor mergers and flyby encounters, 
which suggest that intermediate-to-high mass galaxies 
experience such frequent interactions that they can rarely 
be considered truly isolated (Mailer et al. 2006; Fakhouri 
& Ma 2008; Jogee et al. 2009; Hopkins et al. 2010; Lotz 
et al. 2011; Sinha & HoUey-Bockelmann 2012). There- 
fore, interaction-induced infiows may play a key role in 
molecular gas replenishment for much of the galaxy pop- 
ulation. 

To explore this idea, this paper examines the relation- 
ship between global H2/HI and recent central star forma- 
tion enhancements parametrized by mass-corrected blue- 
ccnteredness (precisely defined in § 2.3.3) for a represen- 
tative sample of galaxies spanning a wide variety of stel- 
lar masses, morphologies, and evolutionary states. We 
find a striking relationship between mass-corrected blue- 
centeredness and H2/HI for nearly all spiral galaxies and 
a fraction of E/SO galaxies, implying a systematic link 
between external perturbations and global H2/HI ratios. 
Intriguingly, our data also reveal that low-mass blue se- 
quence E/SO galaxies - i.e., E/SO galaxies that fall on 
the blue-sequence in color versus stellar mass space with 
spirals (Kannappan et al. 2009) - define an evolutionary 
track offset from the main relation toward stronger blue- 
centered color gradients at a given H2/HI, and trends in 
these galaxies' total gas-to-stellar mass ratios along this 
track suggest the likelihood of fresh gas accretion during 
a late- to post-starburst phase. Thus we find that several 
evolutionary tracks can be summarized within a "fueling 
diagram" that links the immediate fuel available for star 
formation, the total gas, and a metric tracing the recent 
interactions that drive the fueling. 

2. DATA AND METHODS 

This section describes our initial sample drawn from 
the Nearby Field Galaxy Survey and designed to be rep- 
resentative of intermediate mass E-Sbc galaxies, followed 
by the larger literature sample we use to expand our data 
set and confirm our results. We also describe our new 
CO(l-O) and C0(2-l) observations and our methods for 
extracting the useful quantities of gas mass, stellar mass, 
blue-centeredness, and mass-corrected blue-centeredness. 
Finally, we discuss the cuts applied to our sample to en- 
sure robust results. 

All distances are calculated using heliocentric velocity 
corrected to the Local Group frame of reference follow- 
ing the method of Jansen et al. (2000b) and assuming 
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Hq = 70 km s ^ Mpc ^ , except in cases where a more di- 
rect distance indicator was available in NED. 

2.1. Samples 
2.1.1. The Nearby Field Galaxy Survey 

Our primary sample comes from the Nearby Field 
Galaxy Sm-vey (Jansen et al. 2000a, b; see also Jansen 
& Kannappan 2001), a set of ~200 galaxies selected to 
span a broad range of i3-band luminosities and mor- 
phologies. The data products of the original survey in- 
clude UBR surface photometry and optical spectroscopy 
(Jansen et al. 2000a, b; Kannappan & Fabricant 2001; 
Kannappan et al. 2002, 2009). The sample also in- 
cludes extensive supporting data relevant to this study. 
Roughly 90% of galaxies have Sloan Digital Sky Survey 
(SDSS) DR8 optical imaging (Aihara et al. 2011). AU 
have near infrared (NIR) data from the 2 Micron All Sky 
Survey (2MASS; Skrutskie et al. 2006) while 55% have 
deeper Spitzer Infrared Array Camera (IRAC; Fazio et 
al. 2004) imaging, mainly from Sheth et al. (2010), Mof- 
fett et al. (2012) and Kannappan et al. (2013, submitted, 
hereafter K13s). In addition, all galaxies have single dish 
21cm data (Wei et al. 2010a, K13s). 

For this study, we analyze the 35 out of 39 NFGS galax- 
ies with CO data that pass the useability cuts applied to 
our final sample (see § 2.3.7). This subset of the NFGS 
has a stellar mass range of lO*'* — 10^"'^ Af©, morpholo- 
gies ranging from E to Sbc, and diverse states of star 
formation (e.g., starbursting, post-starburst, quiescent). 
Most of the CO data for this sample came from new 
observations, and unlike many previous investigations of 
molecular gas in galaxies, our NFGS CO observations 
were not in general designed to emphasize CO-bright 
galaxies, but were instead designed to reach strong, sci- 
entifically useful upper limits in the case of CO non- 
detections. 

2.1.2. Literature Sample 

To strengthen our results, we supplement our sample 
with galaxies from the literature with available CO, HI, 
and multi-band imaging data. 

Our literature sample includes galaxies from three 
large surveys: the Spitzer Near Infrared Galaxy Sur- 
vey (SINGS; Kennicutt et al. 2003), ATLAS-3D (Young 
et al. 2011), and COLD CASS (Saintonge et al. 2011). 
These surveys are dominated by high mass galaxies, so 
to supplement the low mass end of our data set, we add 
galaxies from Barone et al. (2000), Garland et al. (2005), 
Leroy et al. (2005)^, Taylor et al. (1998), and Kannap- 
pan et al. (2009). Most of these references are themselves 
the sources of the CO data, although some CO data for 
galaxies in SINGS and Kannappan et al. (2009) come 
from Leroy et al. (2009), Albrecht et al. (2007), Sage 
et al. (2007), Zhu et al. (1999), or our own observations 
(§ 2.2). HI data often come from the same source as 
the CO data, or else from alternate sources in the lit- 
erature or HyperLeda (Huchtmeier et al. 1995; Smoker 
et al. 2000; Salzer et al. 2002; Paturel et al. 2003; Gar- 
land et al. 2004; Springob et al. 2005; Meurer et al. 2006; 
Waher et al. 2008; CatineUa et al. 2010; Haynes et al. 
2011; Catinella et al. 2012; Serra et al. 2012). All optical 
data come from the SDSS DR8, except for a subset of 
the SINGS sample outside the SDSS footprint that has 



BVRI photometry (Muiioz-Mateos et al. 2009). NIR 
imaging data is available for all galaxies from 2MASS. 

Combined, these data sources bring our full sample 
(only considering galaxies that have all necessary data) 
to 627 galaxies. However, our total decreases to 323 
galaxies after we institute a number of useability cuts 
(see § 2.3.7). 

2.2. A^ew Molecular Gas Data 

CO(l-O) data, which we use to estimate molecular gas 
masses, already existed for a handful of our NFGS galax- 
ies prior to this study (Sage et al. 1992; Wei et al. 2010b). 
The rest of the NFGS sample was observed with the Insti- 
tut de Radioastronomie Millimetrique (IRAM) 30m and 
Arizona Radio Observatory (ARO) 12m single dish tele- 
scopes. Total integration times were set by how long it 
took to reach reasonable integrated signal-to-noise ratios 
(S/N> 5) or strong upper limits (yielding H2/HI<0.05) 
on the CO flux. Most observations were single pointings 
toward galaxy centers, but offset positions were observed 
for a handful of larger galaxies. 

Initial observing runs on the IRAM 30m took place in 
Fall 2008 and used the ABCD receivers to observe the 
CO(l-O) and C0(2-l) lines at 115 and 230 GHz simulta- 
neously in both polarizations. The 4 MHz filter bank pro- 
vided 1 GHz bandwidth at velocity resolutions of ^^10. 4 
km s^^ and ~5.2 km s~^ at these two frequencies. Fur- 
ther observations were taken in Fall/Winter 2009/2010 
with the newly commissioned EMIR receiver in conjunc- 
tion with the WILMA backend, which supplied 2 MHz 
resolution and a total bandwidth of 3.7 GHz. For all 
observations, wobbler switching was used with a throw 
of 2'and individual scans of 6 minutes. The data were 
calibrated via observations of an ambient temperature 
load. The absolute calibration is accurate to 15-20%. 
The half-power beam widths are 22"and 12"at 115 GHz 
and 230 GHz respectively. 

The ARO 12m observations took place between De- 
cember 2010 and April 2011. We used the ALMA 3mm 
receiver in conjunction with both 2 MHz filterbanks (one 
for each polarization), which provided a total bandwidth 
of 512 MHz. We simultaneously used the Millimeter 
Auto Correlator (MAC), with a resolution of 781.2 kHz 
and a total bandwidth of 800 MHz. Observations were 
carried out in beam switching mode with typical throws 
of 2'-4', and individual scans of 6 minutes. The data were 
calibrated using measurements of a noise diode between 
scans, and galaxies with previous observations were used 
to check the calibration. Some of our ARO time was used 
to obtain CO data for five extra galaxies in our literature 
sample. 

The data were reduced using CLASS^. Scans were aver- 
aged together and any bad channels flagged. The spectra 
were Hanning smoothed to a resolution of 10.4 kms~^. 
Baselines were then fit to emission free regions of the 
spectrum, with polynomials of order <5. 

Integrated fluxes and other measured quantities are re- 
ported in Tables 1 and 2. We convert the IRAM 30m 
data from the measured T^ scale to Janskys using Jy/K 
= 6.12 and 6.3 for the EMIR and ABCD receivers mea- 
sured at 115 GHz, and Jy/K = 7.86 and 7.5 for the 

® We do not include marginal detections due to the authors' 
predictions of a high false positive rate. 
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EMIR and ABCD receivers measured at 230 GHz (sec 
IRAM 30m documentation^). The ARO data are ini- 
tially recorded in the Tp scale, which we then convert to 
Janskys using Jy/K — 38.3 (see the ARO 12m documen- 
tation^). 

Fluxes were determined by summing the channels 
within the line profile. The specific integration ranges 
are given in Tables 1 and 2 and were judged by eye for 
each case. If the profile edges were unclear, we made 
the velocity ranges large enough to ensure all flux was 
included without any doubt and also yield a more con- 
servative error estimate. The uncertainty in the total 
flux measurement is given by 



'^f 



,AV^/N^a 



(1) 



Here, arms is the rms noise of the spectrum in Jy as 
measured from line-free channels, /S.V is the velocity res- 
olution in km s~^, and Nchan is the number of channels in 
the integration. For non-detections, we take upper lim- 
its to be 3(7/ , measured over a velocity range defined by 
the larger of the HI line width or an equivalent linewidth 
from an Ha or stellar rotation curve (see K13s). 

CO linewidths (W,5o) are determined by finding where 
the data are greater than 0.5 times the peak flux minus 
the rms noise for 3 consecutive channels, starting from 
the outside of the emission line and working inwards. The 
final left and right edges for width determination are lin- 
ear interpolations to get the fractional channels where 
the data cross the line height. Heliocentric velocities are 
defined to be midway between the two edges found by 
the above algorithm. Following the examples of Schnei- 
der et al. (1986) and Fouque et al. (1990), we estimate 
line width uncertainties by generating a series of artificial 
observations over a model grid with varying line steep- 
ness and peak signal-to-noise ratios. At our resolution of 
10.4 kms^^, the standard deviation of the linewidth mea- 
surements at each grid point is approximately described 

by 

pO.75 

awKfi — 6-8 TT-Q (2) 

(S/N)°-^ 

where P is the steepness parameter (defined as P = 
{W50 — W20) /2) and S/N is the peak signal-to-noisc ra- 
tio. We stress that linewidths for profiles with peak 
S/N< 6 become extremely unreliable, and should be used 
with caution. We include linewidth measurements in ta- 
bles 1 and 2 for completeness, but these measurements 
never enter into the analysis of this paper. 

2.3. Methods 

2.3.1. Optical and NIR Pfiotometry 

Optical/NIR photometry is needed to estimate stellar 
masses and track recent central star formation enhance- 
ments. For our analysis, we do not use products of the 
SDSS pipeline since it is prone to shredding extended 
sources and docs not make use of the most recent back- 
ground subtraction algorithm (Blanton et al. 2011). We 
instead recalculate total magnitudes using our own cus- 
tom pipeline, described in greater detail by Eckert et al. 

■^ http://www.iram.fr/IRAMFR/GILDAS 

* http://www. iram. es/IRAMES/main Wiki/Iram30mEfficiencies 



(in prep.). After the downloaded images are co-added, 
bright stars and interloping galaxies are masked. The 
masking process is automated with the aid of SExtrac- 
tor (Bertin & Arnouts 1996), but each mask is inspected 
by eye and adjusted when there is clear over- or under- 
masking using the automatic routine. The ELLIPSE task 
in IRAF is used to extract surface brightness profiles of 
constant center, PA, and ellipticity, and to sum up the 
flux within each isophote. For the NFGS sample, we 
adopt the same PAs and ellipticities used for the UBR 
photometry (Jansen et al. 2000b). For our literature 
sample, we adopt the method of Eckert et al. (in prep.), 
who use ELLIPSE to determine the best PA and ellip- 
ticity from the low surface brightness outer disk of each 
galaxy using the coadded gri images. Models of each 
galaxy are created with the resulting surface brightness 
profiles, and are used to fill in masked regions to correct 
the total fiux. 

Total magnitudes are calculated two ways. First, we 
adopt a curve of growth technique very similar to the 
one outlined in Muiioz-Mateos et al. (2009), where the 
outer disk values of the enclosed magnitude and its radial 
gradient are fit with a linear function. The j/-intercept of 
this line (i.e., where the enclosed magnitude is no longer 
increasing) is the total magnitude. Total magnitudes are 
also calculated by fitting the outer disk of the surface 
brightness profile with an exponential function, similar to 
the method in Jansen et al. (2000b), except that outlier 
points are rejected from the fit. The total fiux is summed 
up to the last isophote used in the fit, after which the fit 
itself is used to estimate the remaining outer flux. 

To estimate systematic uncertainties in our total mag- 
nitudes, we perform each of these total magnitude ex- 
trapolations using slightly differently defined fit ranges 
(between 1 and 8 times the sky noise, between 3 and 10 
times the sky noise, within a 1 mag arcsec'^ range end- 
ing at 1 or 3 times the sky noise, and finally using the 
last 5 data points above 1.5 times the sky noise) and then 
average the results (ignoring > 6a outliers) to obtain our 
final magnitude. We take the difference of the maximum 
and minimum magnitude estimates divided by two (also 
ignoring outliers) as our official systematic uncertainty. 
This uncertainty is added in quadrature with the Poisson 
statistical uncertainty. 

The ellipses used for total magnitude calculation best 
match the shape of the far outer disk or halo, and are 
therefore not ideal for calculating galaxy inclinations. To 
estimate photometric inclinations, ELLIPSE is run a sec- 
ond time where ellipticity is allowed to vary while PA is 
kept fixed to the previously used value. The ellipticity 
at a surface brightness of 22.5 mag arcsec"^ in the coad- 
ded gri images usually reliably traces the higher surface 
brightness inner disk and provides accurate inclinations, 
which are estimated using 



cos I 
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Here, q is the intrinsic disk thickness (assumed to be 
q = 0.2), and b/a is the minor-to-major axis ratio de- 
rived from the ellipticity. For the NFGS sample, we use 
the same inclinations as K13s for consistency, which are 

http://aro.as.arizona.edu/12.obs.manual/chapter.3.htm#3..ReceiveffygQf^ O^ ^llC game equation, but not USing the b/a mca- 
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surcd from SDSS photometry. 

JHK magnitudes are also recalculated using the 
2MASS imaging data. Here, we redo the background 
subtraction for the 2MASS images using a method simi- 
lar to that used in the original 2MASS pipeline (Jarrett 
et al. 2000), where the sky was fit by 3rd order polyno- 
mials. However, we fit the background of each relevant 
frame only in the region local to the galaxy of interest 
(within ~ 5 X -R25), with the galaxy itself and any stars 
or background galaxies masked. We impose the param- 
eters from our first set of optical surface brightness pro- 
files (center, PA, and ellipticity) to extract NIR surface 
brightness profiles. 

The two methods for calculating total optical magni- 
tudes described above are again used to calculate the to- 
tal NIR magnitudes. However, tests comparing the out- 
put of these two methods against magnitudes calculated 
from much deeper Spitzer IRAC 3.6/im imaging reveal 
that the exponential fit method performs systematically 
worse than the curve of growth method when applied to 
the relatively shallow 2MASS data. The fits also perform 
the best when using a fit range between 3 and 10 times 
the sky noise. We thus solely use the curve-of-growth 
derived total magnitudes determined using this fit range 
as our final estimates, although we still rely on the dif- 
ference between the curve-of-growth and exponential fit 
derived magnitudes to estimate the systematic error for 
most galaxies. For any galaxy for which our different 
magnitude estimation techniques yield very different re- 
sults (disagreement of more than 0.5 mags), we resort 
to taking an aperture magnitude of the galaxy, and then 
infer the total magnitude by multiplying by the total-to- 
aperture fiux ratio of a higher S/N passband (another 
NIR band like J if possible, otherwise i band). 

All optical and NIR magnitudes are corrected for fore- 
ground extinction using the dust maps of Schlegel et al. 
(1998) and the extinction curve of O'Donnell (1994). 

2.3.2. Stellar Masses 

Stellar masses arc calculated with an improved version 
of the method described by Kannappan et al. (2009). 
Taking as inputs a combination of UBR, ugriz, JHK, 
and 3.6//m photometry, along with global optical spec- 
tra (or any subset of these inputs that are available; 
see K13s and references therein), the stellar mass esti- 
mation code fits mixed young + old stellar populations 
built from pairs of simple stellar population (SSP) mod- 
els from Bruzual & Chariot (2003) assuming a Salpeter 
IMF. Output stellar masses are scaled by 0.7 to match 
a "diet" Salpeter IMF containing fewer low mass stars 
(Bell et al. 2003). The only significant change compared 
to Kannappan et al. (2009) is the addition of a very young 
SSP with age 5 Myr to the suite of models. Note also 
that model SSP pairs can include a "middle-aged" young 
SSP, as long as its age is younger than the old SSP, which 
was also true for the Kannappan et al. (2009) model grid. 
In addition, the UBR zero points have been adjusted for 
consistency (see K13s). The final stellar mass estimate is 
defined by the median and 68% confidence interval of the 
likelihood weighted mass distribution over the full model 
grid. 

Stellar masses for our literature sample are calculated 
using only SDSS and 2MASS photometry. To determine 
whether the lack of spectroscopy leads to any system- 
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Figure 1. Difference between the highest quality stellar masses for 
the NFGS (log M^°=S calculated from UBR, SDSS, 2MASS, and 
IRAC photometry plus integrated spectroscopy) and stellar masses 
for the same galaxies estimated with data available for our litera- 
ture sample (Log M'^'', calculated using only our custom SDSS and 
2MASS photometry), testing the quality of stellar masses for our 
literature sample. The solid horizontal line represents 1:1 agree- 
ment, while the dashed line shows the ordinary least-squares fit 
with outlier rejection, which is not different at a statistically sig- 
nificant level. The scatter has negligible effect on our results. 

atic differences in our stellar mass estimates, we compare 
the high quality stellar masses from the NFGS against a 
second set calculated using only our custom SDSS and 
2MASS photometry. As shown in Figure 1, the two 
methods of estimation are in good agreement. There 
is no statistical significance in the weak linear trend be- 
tween the two mass estimates as function of stellar mass, 
and the la scatter between the two estimates is 0.04 
dex, much less than the typical uncertainty of ~0.2 dex 
for stellar masses. 

2.3.3. Blue-Centeredness 

To track enhancements in recent central star forma- 
tion, we rely on a simple measure of the color gradi- 
ent referred to as blue- center edness (AC), defined as the 
outer disk color from the half-light radius (rso) to the 
75% light radius (^75) minus the color from the center 
to r5o (Jansen et al. 2000b). ^^ For future reference, gen- 
eral discussion of bluc-ccnteredness will use the notation, 
AC, while specific discussion involving a particular color 
will note that color explicitly, e.g. lS.{g — r). The radii in 
this definition reliably separate bulge/disk colors with- 
out being sensitive to variations in bulge-to-disk ratio 
(Kannappan et al. 2004). The half-light radius is a natu- 
ral separator of inner/outer galaxy growth; for example, 
Peletier & Balcells (1996) note shifts in colors and ages at 
approximately this radius. We also stress that this sim- 
ple measure appears robust to variations in the dividing 
radii. For example, using the 40% and 90% enclosed light 
radii returns approximately the same results. 

Before blue-centeredness is used to measure enhance- 
ments in recent central star formation, we account 
for other galaxy properties that may affect the color 
gradients. Kannappan et al. (2004) find that blue- 
centeredness correlates with galaxy luminosity, likely due 

^^ Whenever we have reprocessed SDSS images, blue- 
centeredness is calculated using the second set of ellipse fits, i.e., 
the ones that are used to estimate inclinations, since they best 
track the star-forming inner disks (see § 2.3.1). 
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to the fact that both metalhcity and stellar popula- 
tion age can influence galaxy color gradients, reddening 
galaxy centers with increasing luminosity. These authors 
remove the luminosity trend by subtracting off the fitted 
relation between A{B — R) and luminosity. We follow the 
same approach, performing an ordinary least-squares fit 
(Isobe et al. 1990) of AC against stellar mass (Figure 2). 
We derive the following relations: 



A(u-r) = (-0.114 ±0.025) log M, 
A(ii -g) = (-0.064 ± 0.017) log M^ 
Mg -r) = (-0.049 ± 0.009) log A/* 



(0.983 ±0.249) (4) 

(0.556 ±0.167) (5) 

- (0.417 ±0.095)(6) 



The use of an ordinary least-squares fit of y vs. x is 
crucial, since in order to remove any trend with stellar 
mass, we must minimize the scatter in AC relative to 
stellar mass. Following Kannappan et al. (2004), we cal- 
ibrate this mass correction using the NFGS sample (not 
just our subset with CO data) since it is our most repre- 
sentative data set and has the most robust stellar mass 
estimates. The fit is limited to star-forming disk galaxies, 
identified as galaxies with morphology SO or later and ci- 
ther detected Ha emission extending beyond the nucleus 
or U — K < 4, which helps exclude quenched galaxies 
(Kannappan & Wei 2008, K13s). We reject galaxies with 
known strong AGN and restrict the fit to M* > 10^-^Mq, 
avoiding the low M* tail of the NFGS where all blue- 
centercdness states may not be evenly represented. In 
total, 135 NFGS galaxies were used to derive these rela- 
tions. The residuals of this fit give mass-corrected blue- 
centeredness (denoted by AC™). In simple terms, AC™ 
is a measure of the color difference between the inner and 
outer regions of a galaxy, relative to the typical color dif- 
ference for galaxies at a given stellar mass, where higher 
values imply bluer centers. More physically, this param- 
eter tracks recent central star formation relative to the 
typical star-forming galaxy at a given stellar mass. 

We find that galaxies with and without strong pecu- 
liarities (i.e., signs of recent interactions) as defined by 
Kannappan et al. (2004) have different underlying dis- 
tributions of A{g — rY"- at a 99% confidence level, so 
our definition of mass-corrected blue-centeredness suc- 
cessfully recreates the correlation found by these au- 
thors using luminosity-corrected blue-centeredness (their 
Fig. 8). However, we stress that while our mass correc- 
tion is physically motivated, it is modest, and the distri- 
butions of uncorrected /S.{g — r) for peculiar/non-peculiar 
galaxies differ at confidence level comparable to that for 
distributions using A(g — j-)™. Our general results pre- 
sented in § 3 are still found without the correction in 
place, as discussed in § 3.2.2. 

Within this paper, the total error budget on AC™ in- 
cludes contributions from the measurement error of AC, 
as well as additional uncertainties due to the error in r.go 
and r75, the error in the stellar mass, and the uncertainty 
in the fitted relation between AC and stellar mass. In 
the majority of cases, the Poisson error is dominant. 

We note that optical color gradients have advantages 
over direct star formation rate tracers (e.g., FUV+24//m) 
for the purposes of our study. A practical advantage is 
that the data needed to compute color gradients are read- 
ily available for most galaxies. In addition, the longer 
timcscales over which optical indicators remain sensitive 
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Figure 2. A(g — r) versus log (M,/Mq), illustrating the correla- 
tion between blue-centeredness (AC) and stellar mass. To isolate 
the effect of recent central star formation enhancements, we re- 
move the stellar mass trend using an ordinary least-squares fit, 
minimizing scatter in blue-centeredness relative to stellar mass 
. The residuals of this fit give mass-corrected blue-centeredness 
(AC"), which measures outer minus inner color relative to the 
typical value for a galaxy at a given stellar mass, allowing us to 
trace recent central star formation enhancements above the norm. 



to enhanced star formation are more useful for study- 
ing the extended evolution of galaxies. They enable us 
to analyze the stellar populations beyond the timescale 
of the star formation and gas consumption itself, allow- 
ing us to note the longer term consequences of recent 
star formation episodes. We defer a detailed analysis of 
the timescales of these optical indicators to future work, 
but Figure 11 of Kannappan et al. (2004) shows a sim- 
ple model wherein the blue-centeredness fades on the or- 
der of 0.5-2 Gyr after the central star formation episode 
ends. It should be noted that this model represents a sin- 
gle case, and the evolution of blue-centeredness may vary 
depending on the duration and size of the burst as well 
as the composition of the preexisting stellar populations 
in the bulge and disk. 

For consistency and since the vast majority of our 
galaxies have SDSS data, we always quote AC™ using 
SDSS-equivalent colors. A handful of galaxies from our 
NFGS sample do not have SDSS data, but do have UBR 
photometry. We derive the following conversions be- 
tween the NFGS UBR and ugriz systems using colors 
measured within the B band 25 magarcsec"^ isophote 
for 183 NFGS galaxies: 



g-r = 0.62(5 - 


-i?)-0.18;a= 


= 0.04 


(7) 


u-.g=1.16(C/- 


-B) + im;a = 


= 0.16 


(8) 


u-r = 0.91(C/- 


-i?) + 0.61; cr = 


= 0.21 


(9) 



These relations are useful only for the NFGS due to 
the possibly different UBR zeropoints used in the NFGS 
relative to other photometric studies (see K13s). The 
relations have been corrected for foreground extinction, 
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although as we are measuring color differences, both fore- 
ground extinction corrections and k-corrections cancel. 
Where these conversions are used, their uncertainties are 
incorporated into the final error on AC™. 

2.3.4. Internal Extinction Effects 

Since we make use of optical color gradients, internal 
dust extinction is a concern. A systematic effect from 
dust may manifest itself as a dependence on inclination, 
since higher inclination galaxies show more dust along 
the line of sight. To determine whether internal extinc- 
tion is systematically biasing our results, we run Spear- 
man rank tests on all star-forming disk galaxies in the 
NFGS with measured inclinations >15° after excluding 
AGN. We find no significant correlation between incli- 
nation and AC™. Shown in Figure 3, the same rank 
test for galaxies with stellar masses above ~ lO^'"^ Mq - 
roughly equivalent to the 120 km s~^ rotation velocity 
threshold above which dust lanes become more promi- 
nent in galaxies (Dalcanton ct al. 2004; see also K13s) - 
implies somewhat significant correlations (3%, 5%, and 
4% chance of being random for A(u — r)™, A(m — p)™, 
and A(g — r)™ respectively). Though weak correlations 
may exist, these trends are much smaller than the scat- 
ter (Fig. 3). Thus we conclude that systematic internal 
extinction effects should have a minimal effect on our 
results. 

We should note that while we find only small system- 
atic effects due to inclination, dust within galaxies will 
inevitably alter our color gradients if it is present. There 
are galaxies in our sample for which visual inspection 
has revealed significant dust features. These galaxies and 
their impact on our results are discussed in § 3 and § 4. 

2.3.5. Gas Masses 

To ensure that our data are as uniform as possible, we 
recalculate all gas masses using the same formulae. The 
HI mass is calculated as (Haynes & Giovanelli 1984): 



Afm = 2.36 X l{f'D^F2i 



(10) 



Here, F21 is the measured 21cm fiux in Jy km s~^ and 
D is the distance to the galaxy in Mpc. The molecular 
gas mass is estimated as (Sanders et al. 1991): 



Mh 



1.18 X 10"* 



X 



CO 



3 X 1020 



D'Fco 



(11) 



Here, Fco is the measured CO flux in Jy km s ^. We 
assume a constant Xco of 2 x 10^" cm^^ {y^ ]jjn s~^)^^ 
(Strong & Mattox 1996; Dame et al. 2001); see § 3.1.3 for 
an analysis of how assuming constant Xco may affect our 
results. In further analysis, we multiply all gas masses 
by a factor of 1.4 to account for helium. 

2.3.6. Beam Corrections 

The majority of our CO fiux measurements come from 
single dish telescopes with single pointings at galaxy cen- 
ters. Unlike 21cm beams, which are typically several 
times the size of our galaxies, the CO beams are often 
comparable to or smaller than our galaxies. However, 
since CO scale lengths are typically smaller than optical 
scale lengths, this size mismatch is not always a major 



issue when trying to estimate total H2 mass. By mod- 
eling the beam pattern of the telescope and assuming 
a CO flux distribution, one can estimate the additional 
flux not detected by a single pointing. We apply beam 
corrections to our single-dish CO data following the same 
method as Lisenfeld et al. (2011). The aperture correc- 
tion is defined as the ratio of the total-to-observed CO 
intensity. 



/ — -fcO,total/-fcO,observed 



(12) 



We assume that the CO distribution Iqq follows an 
exponential disk. 



/co(r)-/oe-^/"°° 



(13) 



where /ico is the scale length of the CO emission. The 
total CO intensity is simply the integral of this CO dis- 
tribution to infinity: 



/, 



CO, total 



/■oo 

/ 27r/ore~''/^°°dr = 2t: IqK\ 
Jq 



2 

CO 



(14) 



The observed CO intensity is the integral of the CO 
distribution convolved with the beam pattern. Assuming 
a Gaussian beam with HPBW, 6, this is written: 



/•oo />CXD 

^CO,observed = 4/co / / exp 

Jo Jo 



\Jx^ + ip- 



«co 



X cxp 



ln2 



2y cos i 
6 



dxdy 



A major uncertainty in this calculation is the value 
of ^co- Studies of star- forming late- type galaxies have 
shown hco ^ (0.2 ± 0.05)i?25, where R25 is the B-band 
25 magarcsec~2 isophotal radius (Young et al. 1995; 
Leroy et al. 2008; Schruba et al. 2011). We assume this 
relation holds for all late type galaxies in our sample. 
Whether the same is true for E/SO galaxies is uncertain, 
as their CO scale lengths have not been as thoroughly 
studied in the literature. We examine the CO profiles 
of Wei (2010; see also Wei et al. 2010b), who mapped 
the CO(l-O) distribution in a sample of E/SO galaxies 
using the CARMA array. We find an average value of 
^co = 0.1i?25 with a standard deviation of 0.05i?25. We 
adopt this scale length for all E/SO galaxies in our sam- 
ple. It should be noted that these radial profiles are 
primarily for blue-sequence E/SOs, not the more com- 
mon red-sequence E/SOs, simply because CO is more of- 
ten detected in the former. The two red-sequence E/SO 
galaxies in the Wei (2010) sample have CO scale lengths 
close to 0.05i?25, but two data points are not enough to 
warrant separate definitions of hco for red- and blue- 
sequence E/SOs. 

We assume uncertainties in hco of 0.05i?25 for both 
spiral and E/SO galaxies. The additional uncertainty 
that propagates into the beam-corrected H2 mass de- 
pends on the size of the galaxy relative to the beam. 
Within our final sample (see § 2.3.7) additional uncer- 
tainties are on the order of 5%, though some are as large 
as 25%. We also note that for galaxies in our sample 
for which we have both single dish CO fiuxes and re- 
solved CO maps from Wei (2010), wc adopt the CO 
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Figure 3. Inclination versus AC" for star-forming disk galaxies with i>15° from the NFGS. The full sample shows no correlation between 
these two variables, implying no systematic effect of inclination (and associated dust extinction) on our results. Red crosses mark galaxies 
with stellar masses above 10^'^ Mq, which show a weak correlations between inclination and AC" (3%, 5%, and 4% chance of being 
random for A{n — r), A{u — g), and A{g — r) respectively). The dotted lines show the inverse least-squares fits (i.e., minimizing the scatter 
in the x-direction) to just the more massive galaxies. Any systematic effects due to inclination appear minor. However, irregular dust 
structures, not probed with this test, will inevitably alter color gradients (see § 3.1.1 and § 4.1.1). 



fluxes from the single dish observations and beam-correct 
these fluxes using the measured scale-lengths from the 
CO maps. 

2.3.7. Useahility Cntena 

After identifying from our NFGS-|-literature compila- 
tion all galaxies with the necessary optical/NIR imaging 
as well as 21cm and CO(l-O) flux measurements (627 
galaxies) , we institute a number of usability criteria that 
galaxies must pass before they are included in our fi- 
nal sample. These are: (1) SDSS r-band half-light radii 
larger than 5" to ensure blue-centeredness calculations 
are not compromised by variations in the PSF between 
the different SDSS passbands; (2) minimal beam correc- 
tions to detected CO fluxes, so that the corrected CO 
fluxes are no larger than 1.5 times the measured fluxes 
(equivalent to a change in H2/HI < 0.2 dex), ensuring 
the estimated H2 masses are minimally dependent on the 
assumed model of the CO distribution; (3) strong upper 
limits on total H2 masses, i.e. where CO detections are 
missing but HI detections are not, H2/HI must be < 0.05. 
We do not remove any galaxies with HI upper limits 
from our sample. With these cuts, our NFGS-(-litcrature 
sample totals 323 galaxies. Included in this tally are 
11 galaxies we judged to be highly peculiar/interacting 
whose derived position angles and ellipticities carry little 
meaning, and whose blue-centeredness calculations are 
therefore untrustworthy; we consider their properties in 
§ 3.1.2. We also include an additional 122 "quenched" 
galaxies (two of which are also counted as highly pecu- 
liar), which pass all useability criteria except for having 
upper limits in both HI and H2 mass; these are discussed 
in § 3.2.3. For our NFGS sample, upper limits from the 
IRAM 30m telescope are preferred over the ARO 12m 
telescope since they are always stronger. Weak upper 
limits led to the removal of 4 out of 39 NFGS galaxies 
with CO data from our final sample, though 3 of the 4 
were upper limits from prior CARMA observations (Wei 
et al. 2010b). Where we have both IRAM 30m telescope 



and ARO 12m telescope detections for our NFGS sam- 
ple, we prefer the IRAM 30m data as long as the beam 
correction is less than 1.5. Otherwise, we use the ARO 
12m data. 

It is important to note that while our sample is diverse, 
it is not statistically representative of the galaxy popula- 
tion since our combined data set is subject to whatever 
biases exist in past CO studies. Nonetheless, our use- 
ability criteria, which enforce minimum and maximum 
apparent radii for many galaxies in our sample, do not 
appear to bias us towards galaxies of a certain physical 
size thanks to the wide variety of distances within our 
sample. The distributions of physical sizes of our galax- 
ies before and after we institute our useability criteria 
are not significantly different, as confirmed by a K-S test. 
Our final sample includes some galaxies taken from the 
NFGS, which were in fact originally chosen to be repre- 
sentative of the galaxy population (see § 2.1.1). We will 
use this subset to investigate our results in the context 
of a broadly representative data set. 

Our final sample, including the quenched and pecu- 
liar galaxies, is shown on the u — r color versus stellar 
mass plane in Figure 4. This sample includes some galax- 
ies with M^, < 10^'^ M0, although the blue-centeredness 
mass-correction was calibrated with only galaxies above 
this mass. These lower mass galaxies do not show un- 
usual behavior within our results, and are therefore kept 
in our final sample (see § 3.2.2 for further discussion). 

Our final compiled data set is made available in ma- 
chine readable format in the online edition of our paper. 
A summary of the data is given in Table 3. 

3. RESULTS 

In this section, we describe our findings on the rela- 
tionship between global H2/HI and mass-corrected blue- 
centeredness. These variables are plotted in Figure 5, 
which we hereafter refer to as the "fueling diagram," 
since it links the fraction of gas available as direct fuel 
for star formation and a tracer of recent interactions that 
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Figure 4. Our final NFGS+Utcrature sample plotted in (u — rY 
versus stellar mass space. The dashed line represents the red/blue 
sequence divider, which is based on the analysis of Moffett et 
al. (2013, in prep.), which identifies the midpoint between the two 
sequence peaks at each stellar mass. Our red/blue sequence divider 
is shifted redder than this midpoint by +0.1 mag. The superscript 
i indicates that an internal extinction correction has been applied 
following Moffett et al. 

drive the fueling. 

In § 3.1.1, we examine the basic structure of the fu- 
ehng diagram, which is broken up into three branches. 
The left branch holds most galaxies and shows a pos- 
itive correlation between H2/HI and AC™, consistent 
with the idea that molecular gas content is systematically 
linked to galaxy interactions as traced by mass-corrected 
blue-centeredness. The right and bottom branches are 
well-defined loci that deviate from this expected trend, 
showing first increasing then decreasing mass-corrected 
blue-centeredness as H2/HI decreases. 

In § 3.2, we explore the properties of galaxies within 
each branch. The left branch is largely occupied by a mix 
of barred and unbarred spiral galaxies, with a wide range 
of stellar masses and gas fractions. Conversely, the right 
and bottom branch are almost completely populated by 
low mass E/SO galaxies, specifically blue-sequence E/SO 
galaxies. 

Lastly, in § 3.3 we use differences in optical color gra- 
dients to explore the evolution of galaxies within the fu- 
eling diagram. We find no preferred direction of evo- 
lution along the left branch, while the galaxies on the 
right /bottom branches appear to be evolving in a clock- 
wise fashion back towards the left branch. Intriguingly, 
along this path there are systematic increases in total gas 
content and a transition from E/SO to spiral morphology, 
which may signal a major transformation as these galax- 
ies progress back towards the left branch on the fueling 
diagram. We collect and interpret these findings in § 4. 

3.1. The Distribution of Galaxies in the Fueling 
Diagram 

3.1.1. The Three Branches 

The fueling diagram is shown in Figure 5, plotted using 
mass-corrected blue-centeredness based on m — r, u — g, 
and g — r colors, all of which recreate the same basic 
structure. The regions of parameter space roughly defin- 
ing the left, right, and bottom branches are overlaid on 
the A(5 — r)™ data in Figure 6. 

The majority of our sample falls on the left branch. 



which shows a positive correlation between II2/HI and 
mass-corrected blue-centeredness. There is considerable 
scatter which appears to increase above H2/HI~0.5, pre- 
dominantly biasing galaxies towards lower AC". We 
attribute at least some of this scatter to centrally con- 
centrated dust. Visual inspection has shown several of 
the galaxies that scatter to the left of the main trend 
have distinct dust features (e.g., M82, which lies at 
A(.g - r)™ = -0.14 and H2/HI=3.6). Galaxies in the 
"dusty zone" do not have preferentially high inclination 
(§ 2.3.4). Rather, their centrally concentrated dust may 
reflect a shared evolutionary state (see § 4.1.1). The 
right branch of the fueling diagram, which begins to ap- 
pear at A(g — r)'" ~ 0.05 — 0.1, is far less populated 
but still well defined. It shows a relationship opposite to 
the left branch, where mass-corrected bluc-centeredncss 
increases while H2/III decreases. The bottom branch, 
which encompasses galaxies below Il2/HI~0.06; shows 
no clear correlation between II2/HI and mass-corrected 
blue-centeredness. 

Among the galaxies in our sample with gas detections 
and reliable AC™ measurements, '^75% are on the left 
branch, '-^5% are on the right branch, and ^20% are on 
the bottom branch. However, these percentages are very 
crude since the branches connect with each other and 
their definitions are not exact. We also stress that our 
sample is not statistically representative of the galaxy 
population, largely because of the lack of CO measure- 
ments for low-mass galaxies. Thus, the fractions of galax- 
ies falling on the left, right, and bottom branches in 
Fig. 5 cannot be used to infer the true frequency with 
which galaxies fall on these parts of the fueling diagram. 
If instead we limit ourselves to the more representative 
NFGS subsample, the fractions of galaxies on the left, 
right, and bottom branches are roughly 65%, 10%, and 
25%. But again, this subsample is only approximately 
representative. 

Between these three branches exists a region where 
no data points lie. This hole is clearest in the ver- 
sion of the fueling diagram using A(g — r)™ because the 
right/bottom branches are most evenly populated. For 
this reason, wc default to using the /S.{g — r)™ version of 
plot to show further results. 

3.1.2. Is the hole real? 

The structure seen in the fueling diagram - particularly 
the unoccupied region between the three branches - does 
not appear to be artificially created by our sample restric- 
tions. As described in § 2.1, our full NFGS-|-literature 
sample was unrestricted in terms of stellar mass, gas 
content, and star formation properties. Subsequent re- 
strictions applied were mostly related to the useability of 
the data. The only selection criterion we can reasonably 
relax is our limit on the level of permitted CO beam- 
correction used to estimate flux missed by the telescope 
beam. As discussed in § 2.3.6, we enforce the corrected- 
to-measured CO flux to be less than 1.5. After relaxing 
the corrected-to-measured CO flux to be less than 2, 5, 
and 10, the three branches of the fueling diagram re- 
main well defined (although with increased scatter) and 
no new galaxies appear to invalidate the existence of the 
empty region between the branches. 

It is interesting to examine whether peculiar or actively 
interacting galaxies could potentially fill the empty re- 
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Figure 5. Relationship between H2/HI and mass-correeted blue-eenteredness siiown using tiiree different color gradients. Grey points 
denote lower and upper limits in H2/Hf . To the left of the y-axis are dots representing H2/HI values for galaxies with pceuliar morphologies, 
many clearly interacting. Due to their disturbed state, we cannot measure AC™. Below each x-axis are dots showing the measured AC™ 
for quenched galaxies. All three panels show distinct loci and an empty region between them. 

3.1.3. Xco Effects 

The structure of galaxies in the fueling diagram ap- 
pears robust against possible variations of Xco due to 
its dependence on metallicity. To test this, we use 0/H 
to estimate Xco for each galaxy, employing the calibra- 
tion from Obreschkow & Rawlings (2009). Estimates of 
0/H from optical line ratios are only available for our 
NFGS (Kewley et al. 2005) and SINGS (Moustakas et al. 
2010) samples and the studies of Barone et al. (2000) and 
Taylor et al. (1998). To be able to carry out this anal- 
ysis for our full sample, we also estimate X^q for each 
galaxy from B band luminosity, again using the calibra- 
tion of Obreschkow & Rawlings (2009), who find it to be 
the next most reliable estimator of Xco after metallicity. 
The values of 12+logO/H in our sample (where known) 
range from 7.9-9.2, yielding Xco estimates between 0.4 
and 10 times the Milky Way value. The calibration using 
B-band luminosity yields a similar range of Xco • 

Figure 7 shows the effect that variable Xco due to 
metallicity has on the fueling diagram. Even with the 
variation in Xco j ths right/bottom branches remain dis- 
tinct from the left branch and the hole remains intact, 
although some of the upper limits on the bottom branch 
show significant increases in H2/HI due to their estimates 
of Xco being several times the Milky Way value. These 
galaxies are mainly dwarfs with stellar masses < a few 
xlO^M0. where this sort of deviation might be expected. 
In general however, the structure of the fueling diagram 
remains unchanged. 

The above analysis ignores other physics that may al- 
ter Xco- Most relevant to our study is the possible de- 
crease in Xco in very high gas surface density regimes 
where the the ISM turns almost entirely molecular. Such 
a situation could occur as the result of infiows that drive 
large amounts of gas to the centers of galaxies, and has 
been observed in many systems from ULlRGs to ordi- 
nary spirals (Downes et al. 1993; Regan 2000). The gas 
in these extremely high surface density regions is ex- 
pected to have increased temperature relative to molec- 
ular clouds in less dense environments (Narayanan et al. 
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Figured. Empirically defined branches in H2 /HI vs. AC" space. 
These data points are the same as in Figure 5c. 



gion. These galaxies lack reliable AC measurements, so 
we do not plot them within the fueling diagram. How- 
ever, we mark their measured H2/HI to the left of the 
y-axis in Figure 5. Their H2/HI values do not exclusively 
cluster in the range spanned by the hole, although some 
have the proper values to fall within it. We conclude that 
if galaxies like those represented in our sample ever fall 
in the empty region, they must do so only briefiy dur- 
ing a phase of rapid evolution. The specific reasons why 
galaxies rarely settle in this region of parameter space are 
not immediately apparent. Detailed modeling of galaxies 
may be necessary to explain this phenomenon and will 
be the focus of future work. 
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CO(2- 



Fr 



CO(2- 



-l)/^CO(l- 

-i)/-Pco(i- 



Figure 7. The possible effect of variations in Xco oil our data, with Xco estimated separately for each galaxy using _B-band luminosity 
(left) and log(0/H) derived from nebular emission lines (right). Lighter arrows denote galaxies with H2/HI upper or lower limits, and X's 
denote galaxies where log(0/H) is unavailable. Some upper limits show signifieant adjustments, but Xco variations do not ehange the 
overall appearanee of three distinet branches. 

this division is useful since the presence or absence of ex- 
tended spiral arms represents a basic transition in struc- 
ture likely strongly linked to star formation history. 

Wc find a clear bimodality in the distribution of E/SO 
and spiral morphologies. Spiral galaxies almost exclu- 
sively fall on the left branch, and in fact, a Spearman 
rank test on all spiral galaxies with H2/HI>0.06 (to avoid 
confusion with the bottom branch) confirms a correla- 
tion between H2/HI and AC™, with ~5cr confidence (for 
u — r, u — g, and g — r, the probabilities of the distribu- 
tions being random are 5 x 10~*, 8 x 10"'', and 1 x 10~^ 
respectively). E/SOs show a much more diverse distri- 
bution throughout the fueling diagram: some occupy 
the left branch with spirals, but also and more notice- 
ably the rest almost completely define the right branch 
and much of the bottom branch. Several of the E/SOs 
on the right/bottom branches are also centrally concen- 
trated Blue Compact Dwarf galaxies (BCDs), which are 
lumped with traditional E/SOs in our simplified morpho- 
logical classification scheme. Towards the left end of the 
bottom branch, there is a strong shift in morphology from 
E/SO to spiral. 

Barred galaxies are also shown in Figure 8a. Except 
for galaxies that have existing bar classifications from 
the NFGS or Nair & Abraham (2010), bar classifica- 
tions were done by eye using SDSS i-band imaging, since 
bars are best identified in bands that trace the stellar 
light (Eskridge et al. 2000). Bars are noticeably ab- 
sent from all galaxies above A(g — r)™ ^ 0.15, which 
includes most of the right branch, as well as sections of 
the left and bottom branches. Within the left branch 
alone, a Spearman Rank test does not suggest a smooth 
correlation of bar fraction with AC" or with H2/HI. Al- 
though our full sample is not statistically representative 
of the galaxy population, we speculate that the same 
processes that produce extreme blue-centcrcdness may 
destroy bars, while milder processes associated with evo- 
lution along the left branch do not. If we limit our ex- 
amination of bars to only the more representative NFGS 



with mass-corrected blue-centeredness and H2/HI, but 
we find no correlations. Any evidence of increased cen- 
tral gas temperature is likely being washed out in our 
global CO measurements. Even if some of the increase 
in H2/HI with AC™ along the left branch is the result 
of overestimated Xcoi this is still consistent with our 
physical interpretation that the left branch of the fueling 
diagram is the result of galaxy interactions and inflows 
(see § 4.1.1). 

For the remainder of this paper, we use molecular gas 
masses estimated as described in § 2.3.5, assuming con- 
stant Xco- 

3.2. Distribution of Galaxy Properties Within the 
Fueling Diagram 

Having established the basic structure of the fueling 
diagram and its reliability, we now explore how galaxy 
properties - specifically morphology, the presence of a 
bar, stellar mass, blue versus red sequence, and gas con- 
tent - distribute themselves throughout the fueling di- 
agram. Doing so may provide an understanding of the 
physical processes that drive the observed trends. The 
galaxy properties discussed in this section are overlaid on 
the fueling diagram in Figure 8 and briefly summarized 
in § 3.2.6. 

3.2.1. Morphology 

Figure 8a displays galaxy morphologies within the fuel- 
ing diagram. All galaxies are classifled by eye using SDSS 
g-band images, but we check our classifications against 
previously published types when available. Galaxies are 
separated into two categories: E/SOs (including SOa) and 
spirals (Sa-Sd). Using the distinction between SOa and 
Sa as the separation between early- and late-type galax- 
ies may be somewhat sensitive to classification error, but 
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Figure 8. (a) Distribution of early (E/SO - SOa) and late (Sa-Sd) morphologies and bars. (6) Distribution of stellar masses in three 
fundamental mass regimes: stellar mass above Mi,, between Mt and Afj, and below Mt, where M;, is the bimodality mass (IO^^'^Mq) 
and Mt is the gas-richness threshold mass (10^'^Mq). Also marked are galaxies that fall below the range of our blue-centeredness mass- 
correction. See § 3.2.2 for more details on these mass regimes. Galaxies with stellar masses outside the mass regime used to calibrate the 
blue-centeredness mass correction are also noted (c) Distribution of red and blue sequence galaxies (d) Distribution of Mhi+H2-I-Ho/A^*. 

forming with disk-like morphology to typically non-star- 
forming with spheroidal morphology (Kauffmann et al. 
2003), (2) M* < Mt = 5 X IO^Mq, where Mt is the 
gas-richness threshold mass, below which there is a sig- 
nificant increase in gas-dominated galaxies (Kannappan 
et al. 2009, K13s), and (3) Mt < M, < Mb, where 
bulged spirals with intermediate mass content arc the 
norm (K13s). Galaxies with A/* < 10*-^ Mq are also de- 
noted with an extra "+" symbol in this figure. These 
galaxies technically fall outside the stellar mass range 
used to calibrate the blue-centeredness mass-correction 
(§ 2.3.3), and Kannappan et al. (2004) suggest their pos- 
sibly irregular patterns of star-formation may make color 
gradients hard to interpret in relation to higher mass 



subsample, we still see no significant trend between AC™ 
and bar fraction within the left branch. Bars are impor- 
tant galactic structures to put into the context of our 
study, since like galaxy interactions, they are thought 
to enable infiows of gas to the centers of galaxies. We 
discuss bars and interpret their role in § 4.3. 

3.2.2. Stellar Mass 

Figure 8b displays the distribution of stellar masses 
within the fueling diagram. Instead of examining 
the continuous distribution of stellar masses, we di- 
vide the data into three characteristic mass regimes: 
(1) A4 > A/b = 3 X IQIOMq, where Mb is the bi- 
modality mass, a stellar mass scale above which the 
population of galaxies goes from being typically star- 
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galaxies. However, we include them in our plots because 
they show no unusual distribution relative to the rest of 
the sample below Mf. Among these galaxies, the most 
extreme blue-centcrcdncss is found for NGC 3738, with 
A(g— r)™ — 0.19. This galaxy is a centrally concentrated 
BCD, making it completely consistent with its neighbor- 
ing galaxies in the fueling diagram. 

The different stellar mass regimes display patterns 
within the fueling diagram. While galaxies in all three 
regimes span the full range of H2/HI, we find almost no 
galaxies above Mf, on the right/bottom branches. In- 
stead, most of the galaxies on these branches fall below 
Mf. There is also a tendency for galaxies below Mt to 
cluster in the lower-left corner of the fueling diagram 
(many having H2/HI upper limits), but elsewhere along 
the rising branch they appear to spread roughly evenly, 
as do galaxies in the higher stellar mass regimes. We 
note that with our low number statistics, the apparent 
tendency of galaxies on the upper right branch to have 
stellar masses between Mt and M^ is not significant and 
should not be over-interpreted given that the range be- 
tween Mt and Mf, is very narrow, and stellar mass esti- 
mation involves typical errors of ^^0.2 dex. 

In § 2.3.3, we described our motivation and meth- 
ods for subtracting the underlying dependence of bluc- 
centeredness on stellar mass, yielding mass-corrected 
blue-centercdness, or AC". To determine how depen- 
dent our results are on this mass correction, we plot the 
fueling diagram using the uncorrected A(g — r), rather 
than A{g — r)™, in Figure 9. The general structure of the 
fueling diagram remains intact, specifically the presence 
of three branches with a hole between them. The distri- 
bution of different mass regimes illustrates the tendency 
of high-mass galaxies to have more red-centered color 
gradients, which originally motivated our mass correc- 
tion. 

3.2.3. Red and Blue Sequences 

Figure 8c shows the distribution of red and blue- 
sequence E/SOs and spirals in the fueling diagram, which 
are classified based on their positions within the u — r 
versus stellar mass plane shown in Figure 4. The left 
branch is composed of a mix of red and blue sequence 
galaxies. Intriguingly, the right and bottom branches are 
almost completely dominated by blue sequence galaxies, 
despite the common assumption that E/SO galaxies al- 
ways fall on the red sequence. High-mass, red-sequence 
E/SO galaxies in our sample are often quenched, i.e., have 
no detected CO or HI emission and total gas-to-stellar 
mass ratio upper limits less than 0.04Mq. These systems 
are shown below the a;-axis in Figures 5, 8, 9, and 10.^^ 

3.2.4. Gas Content 

Figure 8d shows the distribution of total gas-to-stellar 
mass ratio (total gas = HI+H2 with a 1.4x mass correc- 
tion for He) within the fueling diagram. Broadly speak- 
ing, Hl-to-stcUar mass ratios and total gas-to-stcUar 

^^ We note that the quenched population does not center around 
AC™ ~ 0, but rather around A{g — r)™ ~ 0.05 (as one color 
example). The A{g — r)™ values are not due to the presence of 
excess recent central star formation, but rather the fact that the 
mass correction applied to the color gradients was calibrated on 
star-forming disk galaxies, which tend to show more red-centered 
color gradients than do passively evolving galaxies. 
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Figure 9. The fueling diagram plotted using blue-centeredness 
without a mass correction. Without the mass correction applied 
to blue-centeredness, there is a tendency for higher stellar mass 
galaxies to have more red-centered color gradients (motivating the 
mass correction), but the underlying structure of the fueling dia- 
gram remains intact. 

mass ratios on the left branch are both anti-correlated 
with H2/HI. However, there is significant diversity be- 
low H2/HI'^0.5, where galaxies with both very high and 
very low gas fractions fall. Above H2/HI'^0.5, the gas 
fractions on the left branch are consistently low. 

Most of the galaxies comprising the right and bottom 
branches have substantial gas content, between 10-100% 
of their stellar mass. Along the bottom branch, both HI- 
to-stellar and total gas-to-stellar mass ratios appear anti- 
correlated with mass-corrected blue-centeredness, that 
is, gas content is higher on the left, more red-centered 
side. Using a Spearman rank test on all unqucnchcd 
bottom branch galaxies (defined as H2/HI < 0.06), we 
find the negative correlation of total gas-to-stellar mass 
ratio and A(g — r)™ has a 1.6% chance of being ran- 
dom. This probability drops to 0.2% when restricting 
the test to M* > 10^'^ Mq (i.e., the mass above which 
the blue-centeredness mass correction was calibrated) . A 
few red-sequence galaxies with very low gas fractions also 
lie on the bottom branch, which are among the minor- 
ity of massive (A/* > Mt) E/SOs whose gas data are not 
upper limits. 

3.2.5. AGN 

While our sample lacks uniform/complete data for 
AGN classification, we display known AGN in Figure 10 
and briefly discuss them for two reasons. First, it is im- 
portant to note that AGN themselves are not the cause 
of strongly blue-centered color gradients, as their light 
contribution is typically too small to have any signifi- 
cant effect on mass-corrected blue-centeredness, which 
uses colors measured over large regions of galaxies. Sec- 
ond, we note that AGN are predominantly seen among 
the high mass galaxies in our sample (sec also Fig. 8b), 
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Figure 10. Known AGN (stars) within tlio fueling diagram. This 
classification is not complete or uniform in our sample due to the 
limited number of classifications available in the literature and in- 
complete nuclear spectroscopy. 

with 25 out of 27 AGN hosted by galaxies with M* > M^. 
Their presence, particularly among the high mass ellipti- 
cal galaxies which show AGN in the quenched regime but 
also among galaxies with detected gas, affects our physi- 
cal interpretation of the fueling diagram, and is discussed 
in § 4.1.1 and § 4.1.3. 

3.2.6. Properties Summary 

To summarize, the main properties of each branch of 
the fueling diagram are as follows: 

The left branch is where spiral galaxies are primar- 
ily found. Many galaxies here are barred, but without 
any discernible pattern with respect to AC™ or H2/HI. 
Galaxies on the left branch cover the full range of stellar 
masses, along with a wide range of HI- and total gas- 
to-stellar mass ratios, which broadly behave inversely to 
H2/HI. 

The right and bottom branches are almost entirely 
populated by unbarred E/SO galaxies, although this type 
distribution transitions into largely spirals on the left side 
of the bottom branch, which is also where we begin to 
see barred galaxies again. Stellar masses fall predom- 
inantly below the bimodality mass, and most fall be- 
low the threshold mass. Gas fractions are moderate to 
high, between 10% and 100% of the stellar mass, and 
show a negative correlation with mass-corrected blue- 
centeredness along the bottom branch. 

We also note that the right /bottom branches are domi- 
nated by blue-sequence E/SO galaxies, i.e., galaxies with 
spheroidal morphologies that fall on the blue sequence 
in color versus stellar mass space, which are thought to 
represent a transitional phase (Kannappan et al. 2009). 
This fact, as well as the properties of galaxies within the 
different branches of the fueling diagram, may reveal im- 
portant clues as to how these branches relate to different 
evolutionary states, as discussed in § 4. 



3.3. Do Galaxies Evolve Within the Fueling Diagram? 

If mass-corrected blue-centeredness and II2/III are 
time-varying properties, galaxies should evolve within 
the fueling diagram. To determine whether there is any 
unified direction of evolution within the fueling diagram 
or any one of its branches, we search for systematic off- 
sets between A{u — g)™- and A{g — r)™. This comparison 
is useful because u—g color reddens faster than g—r color 
after a star forming event, and therefore A(m— g)™ tracks 
recently enhanced central star formation over shorter 
time scales and will return to low values faster than 
A(g — r)™. Before comparing these two measurements, 
we must first account for the fact that the range of val- 
ues for A(u — g)™ is larger than for A{g — r)™ (see Fig- 
ure 5). We therefore divide A(u — g)™ and A{g — r)"^ 
by their median absolute deviations (0.102 and 0.054 re- 
spectively, found with the same NFGS subsample used 
to derive Eqs. 4-6) to obtain normalized mass-corrected 
blue-centeredness values, A{u — g)"^* and A(5 — ?')'"*, for 
which the values cover a similar range while not chang- 
ing the position of zero (i.e., the average gradient for the 
population) . 

In Figure 11, arrows display the relative positions of 
A{u - g)™* (head of arrow) and A{g — r)™* (tail of ar- 
row). Galaxies on the left branch are directed randomly 
left or right along it. This result suggests that galax- 
ies either do not evolve throughout this region of the 
fueling diagram, or the evolution is not necessarily in a 
unified direction. Given the prior association of mass- 
corrected blue-centeredness with interactions (Kannap- 
pan et al. 2004; Gonzalez-Perez et al. 2011), we suppose 
that galaxies may oscillate along the left branch, rising 
when infiows enhance central star formation and II2/HI 
(with some of the apparent rise in H2/HI possibly as- 
sociated with a decrease in Xco due to the central gas 
concentration), and fall along the same locus as outer 
disk gas and star formation are renewed. 

Conversely, the loop composed of the right and bottom 
branches shows some partial systematic behavior. Al- 
though the arrows on the right branch show no preferred 
direction, the bottom branch shows an excess of galaxies 
with A(u — 5)™* lower than A{g — r)™*, implying that 
these galaxies are evolving leftwards towards lower mass- 
corrected blue-centeredness. Ignoring the bottom-left re- 
gion of the fueling diagram where the left and bottom 
branches cannot be distinguished (A(g — r)™ < 0.05), 
we find the frequency of leftward facing arrows on the 
bottom branch is higher than the frequency found on 
the left branch at the 4.2(T confidence level, assuming the 
uncertainty in the number of rightward/leftward facing 
arrows follows Poisson statistics. The right and bottom 
branches are likely closely linked (discussed further in 
§ 4), leading to our physical interpretation that galaxies 
on the right branch are still undergoing a central star- 
burst, after which they move to the bottom branch where 
their young central stellar population ages and fades. If 
star formation is still progressing progressing on the right 
branch, then there may be no reliably predictable differ- 
ence between A{u — g)"^* and A{g — r)"**, which would 
explain the lack of a unified direction of arrows for right 
branch galaxies in Figure 11. At this point, we do not 
have an estimate of the timescales associated with the 
evolution along the the right or bottom branches, but 
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Figure 11. Arrows going from A{g — r)'"* to A{u — 3)'"*. Since 
u~ g color is more sensitive to high mass short Uved stars, central 
enhancements in u — g color should be shorter lived than central 
enhancements in 3 — r color. The data at the bottom branch sup- 
port this expectation by showring A(g — r)""* more blue-centered 
than A{u — g)"^* in most cases, indicating that these galaxies' cen- 
tral star formation has ceased and the young population is fading. 
This color comparison suggests a uniform evolutionary direction of 
galaxies within this region of parameter space. 



estimating these timescales by comparing the galaxy col- 
ors to stellar population synthesis models will be a focus 
of future work. We also note that the results of Fig- 
ure 11 are not dependent on the stellar mass correction 
applied to blue-centeredness as the same result is found 
even without any mass correction whatsoever. 

Since comparing A{g — r)™* and A(m — g)™* gives a 
clear direction of evolution along the bottom branch, we 
can tell that galaxies here shift from primarily E/SO to 
primarily spiral morphologies, as well as generally toward 
increasing gas-to-stcllar mass ratios. We interpret this 
pattern as a sign of disk rebuilding in § 4.2. 

4. DISCUSSION 

In § 3 we presented the fueling diagram relating global 
H2/HI to recent central star formation enhancements. 
Having described the three-branch distribution of galax- 
ies within the fueling diagram, the variation of galaxy 
properties along and between the branches, and the ap- 
parent evolution within the diagram, we now collect these 
results to provide an interpretation of the physical pro- 
cesses that drive it. We first discuss the role of galaxy 
interactions in driving H2/HI ratios, considering the im- 
portance of the merger mass ratio, stellar mass, and 
gas richness of the galaxies involved. We then explore 
how trends within the fueling diagram support a sce- 
nario of fresh gas accretion and stellar disk rebuilding 
along the bottom branch. We finish by reassessing the 
validity of the assumed link between mass-corrected blue- 
centeredness and galaxy interactions in light of our new 
results, specifically addressing the role of bars. 



4.1. Mergers and Interactions in the Fueling Diagram 

In the following section, we describe the role that merg- 
ers play in driving the evolution in each branch of the 
fueling diagram. 

4.1.1. The Left Branch and Regions Above It 

The very existence of the left branch provides support 
for the idea that local galaxy interactions play a key role 
in replenishment of molecular gas: since mass-corrected 
blue-centeredness is a signpost of recent galaxy interac- 
tions (Kannappan et al. 2004), and it is correlated with 
H2/HI, galaxy interactions appear to be linked to H2/HI 
in a systematic way. The analysis of A{u — g)™* and 
A(g — 7-)™* in § 3,3 implies that galaxies on the left 
branch may evolve along it in both directions: up as 
a result of each galaxy encounter, boosting the galaxy 
to higher H2/HI and AC™ (with a possible contribution 
due to overluminous CO), then down as the molecular 
reservoir is consumed and the young central stellar pop- 
ulation fades. Total gas-to-stellar mass ratios support 
this scenario. We expect HI and total gas-to-stellar mass 
ratios to decrease at higher AC" since the gas is be- 
ing converted into H2 and then stars, but we also expect 
a mix of gas fractions at low AC" since galaxies settle 
here before or after each burst event while accreting fresh 
disk gas. The bulk of the evolution along the left branch 
must be driven by minor rather than major mergers or 
interactions, since most galaxies appear to retain their 
spiral morphologies. An alternate possibility is that the 
relation along the left branch reflects bar-driven inflows. 
However, bars do not show any statistically significant 
correlation with AC™ along the left branch. The possi- 
ble role of bars is explored further in § 4.3 

Major mergers in the high stellar mass and/or low 
gas fraction regime may help to explain the E/SOs at 
high H2/HI on the left branch (in contrast, gas-rich, low 
mass, but comparable mass ratio mergers appear to fol- 
low a different path; see § 4.1.2-4.1.3). Alternatively, 
these galaxies could be formed by repeated minor merg- 
ers rather than a single event (Bournaud et al. 2007), 
and they could therefore be galaxies that have made 
several oscillations along the left branch. Either way, 
many of the E/SOs at high H2/HI on the left branch are 
likely moving into the quenched regime, up and off the 
plot, as their remaining gas reservoirs have become al- 
most entirely molecular and may soon be completely de- 
pleted. In addition, some of the high stellar mass E/SOs 
at the peak of the rising branch may have previously 
been quenched but recently experienced small gas ac- 
cretion events, possibly associated with satellite accre- 
tion (e.g.. Martini et al. 2013). Some of these E/SOs 
host AGN, as do many of the high mass E/SOs in the 
quenched regime. A gas accretion event could provide 
fuel for AGN activity coupled with a central molecular 
gas concentration, resulting in a high H2/HI ratio with 
minimal change to AC™ (since most of the star forma- 
tion would be occurring very close to the AGN). Galaxies 
undergoing this process would jump vertically in the fuel- 
ing diagram, between the quenched regime and the peak 
of the left branch. This path could place them in the 
hole of the fueling diagram, although they would likely 
move through it relatively quickly. 

As previously noted in § 3.1.1, there is substantial scat- 
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tcr above the left branch towards more red-centcrcd color 
gradients, hkely caused by internal extinction. While 
dust effects may be altering the measured color gradi- 
ents, they usefully highlight galaxies in early stages of 
star-forming events when the young stars are still heav- 
ily embedded in dust clouds. We suspect galaxies in this 
"dusty zone" represent a stage very soon after mergers or 
interactions that are mild enough not to have driven the 
galaxy into the peculiar morphology region off the plot, 
shown on the left hand side of the fueling diagram. As the 
star formation progresses, the dust is likely to clear, al- 
lowing these galaxies to develop the bluer-centered color 
gradients expected for their H2/HI ratios. 

4.1.2. The Right Branch 

We suggest galaxies on the right branch may be the 
result of gas-rich mergers, specifically between galaxies 
of comparable stellar mass. Along with their high val- 
ues of mass-corrected blue-centeredness (suggestive of a 
strong central starburst), the spheroidal morphologies of 
galaxies on the right branch arc consistent with recently 
violent histories. Most of these galaxies are classified as 
blue-sequence E/SOs, and several are also classified as 
BCDs (e.g., Haro 2, NGC 7077, UM 465). The existence 
of blue sequence E/SOs and BCDs in the same regime of 
the fueling diagram argues in favor of them experiencing 
similar evolutionary processes. Previous observational 
and theoretical studies of BCDs (e.g., Ostlin et al. 2001; 
Pustilnik et al. 2001; Bekki 2008) also support merger 
driven evolution. 

While most of these galaxies do not show obvious out- 
ward signs of a recent strong interaction in their opti- 
cal images, such as irregular structure or tidal features 
(the lack of these features is actually built in to our 
analysis since we do not plot highly peculiar/interacting 
systems), smooth optical morphology is not inconsistent 
with a merger having recently occurred. Merger simu- 
lations find the strongest morphological disturbances be- 
fore the peak of induced star formation, which in turn 
typically occurs before galaxies land on the right branch, 
and the complete coalescence of the two merging galaxy 
nuclei commonly happens several hundred Myr before 
the main starburst event (Lotz et al. 2008, 2010). The 
diverse arrow directions in Figure 11 on the right branch 
suggest these galaxies are still actively forming stars well 
after the merger remnant has settled. Signatures of re- 
cent mergers may be more obvious in observations of gas 
morphology and kinematics. HI maps can be extremely 
useful since they trace extended structure and can re- 
tain signatures of interactions as long as 1 Gyr after the 
events (Holwerda et al. 2011). For example, the high 
resolution HI map of Haro 2 (Bravo- Alfaro et al. 2004) 
shows the HI kinematic and optical major axes to be al- 
most perpendicular, consistent with a merger or recent 
accretion event. 

Blue-sequence E/SOs are known to emerge primarily 
below Mb, and become abundant below Mt (Kannappan 
et al. 2009). As seen in Fig. 8, the blue E/SOs on the 
right branch are consistent with this pattern. However, 
the existence of the right branch cannot be driven by 
stellar mass alone. Low stellar mass galaxies are found 
throughout the fueling diagram, and if AC"" were dic- 
tated solely by stellar mass (i.e., equal size bursts occur- 
ring in higher/low mass galaxies yielding lower/higher 



AC™), then the hole seen in the fueling diagram should 
not exist. A merger origin for the right branch is more 
consistent with such a large gap between the left and 
right branches. Furthermore, gas richness likely produces 
distinct evolutionary tracks for galaxies: gas-rich merg- 
ers drive galaxies along the right branch, while gas-poor 
mergers drive galaxies into the quenched regime, up and 
off the plot, as discussed in § 4.1.1. The association of the 
gas-rich merger track with galaxies below Mf, and espe- 
cially Mt is a simple consequence of increasing gas rich- 
ness below those scales (Kannappan et al. 2009, K13s). 

4.1.3. The Bottom Branch 

The bottom branch appears to be part of the same 
evolutionary sequence as the right branch but at a later 
stage. Galaxies here show many of the same properties 
as those on the right branch in terms of stellar mass, gas 
richness, and prominence of blue-sequence E/SOs. How- 
ever, they show depressed H2/HI ratios while uniformly 
evolving leftwards in the fueling diagram (§ 3.3). These 
observations are all consistent with a scenario where 
these are post-starburst galaxies with depleted central 
gas and fading young central stellar populations. 

4.2. Evidence for Disk Rebuilding 

We have argued that relatively low mass, gas-rich, but 
roughly equal-mass-ratio mergers appear to be responsi- 
ble for the creation of blue-sequence E/SO galaxies on the 
right and bottom branches of the fueling diagram. Com- 
bining the known direction of evolution on the bottom 
branch (§ 3.3) with the observed trends in gas-to-stellar 
mass ratio and morphology (§ 3.2.1,3.2.4) further implies 
that these galaxies may regrow gas and later stellar disks. 
For galaxies on the bottom branch, there is a general 
increase in the total gas-to-stellar mass ratio as AC™ 
decreases. Since galaxies here appear to be evolving left- 
ward on the fueling diagram, their total gas content must 
be growing as their central stellar populations fade. In 
the same direction, there is a transition from E/SO mor- 
phologies to spiral morphologies. These combined trends 
arc consistent with a scenario of fresh outer-disk gas ac- 
cretion and eventual conversion into visible spiral arms, 
and in fact blue-sequence E/SOs have ideal stellar surface 
mass densities for turning gas efficiently into stars, pro- 
moting stellar disk rebuilding (Kannappan et al. 2009; 
see also Kauffmann et al. 2006). Notwithstanding this 
self consistent picture of morphological transformation, 
there remain a handful of blue sequence E/SOs in the bot- 
tom left corner. One possible explanation for their pres- 
ence is that spiral structure formation has been inhibited. 
Two examples where such inhibition may be occurring 
arc described in Kannappan et al. (2009): NGC 7360 
(A(.g - r)'"=0.027, H2/HRO.O23), which hosts counter- 
rotating stellar disks, and UGC 9562 (A(5-r)'" =-0.002, 
H2/HI=0.011), which is a polar ring galaxy. In both of 
these cases, peculiar kinematics may be stifiing spiral 
arm formation. 

Falling below the gas-richness threshold mass Mf 
(§ 3.2.2, K13s) may enable galaxies to evolve leftward 
along the bottom branch. Bottom branch galaxies typ- 
ically fall below M^ and have not only high gas frac- 
tions (as is typical for blue-sequence E/SOs in general, 
Kannappan et al. 2009, Wei et al. 2010a), but gas frac- 
tions that increase as their central starbursts fade. The 
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high/increasing gas fractions may indicate rapid gas ac- 
cretion below Mt as argued by K13s. Theoretical stud- 
ies of gas accretion from the last decade might suggest 
that M( reflects the critical mass scale for cold-mode 
accretion (Birnboim & Dekel 2003; Keres et al. 2005). 
However. Nelson ct al. (2013) calls this interpretation 
into question, arguing that there is no strong transition 
from cold to hot mode accretion, and that the amount 
of gas accreted via the cold mode is not as large as pre- 
viously thought. Regardless of the mode of accretion, 
the high/increasing gas fractions on the bottom branch 
may be explained using the halo mass dependence of gas 
cooling times (Lu et al. 2011). 

Conversely, galaxies on the bottom branch that have 
abnormally low gas fractions relative to the rest of the 
population have masses above M^. One possible rea- 
son may be reduced accretion. Above the bimodality 
mass M;, in particular, both observations and theory sug- 
gest significantly quenched cosmic accretion onto galaxies 
(e.g., Gabor & Dave 2012, K13s). Another possible rea- 
son certain galaxies might fail to accrete fresh gas on the 
bottom branch is their environments. Galaxies in dense 
clusters and groups have long been observed to have de- 
pressed HI fractions (Giovanelli & Haynes 1985; Solanes 
et al. 2001; Cortese et al. 2011), and galaxies near massive 
companions are also more likely to be quenched, even in 
the dwarf regime below Mt (Geha et al. 2012). One likely 
example of neighbor-inhibited accretion within our sam- 
ple is NGC 3073, a blue-sequence E/SO below Mj that 
has an abnormally low gas-to-stellar mass ratio, but lies 
very close to its much larger companion, NGC 3079. The 
fact that NGC 3079 has an AGN and observed outflows 
(Cecil et al. 2001) may also be related to the low gas frac- 
tion in NGC 3073. Other outliers on the bottom branch, 
NGC 4111 and NGC 4270, reside in dense groups. A full 
environmental analysis has not been performed on our 
sample, but these anecdotal cases hint that environment, 
as well as stellar mass, likely determines which galaxies 
are capable of regrowing disks. 

A final possible explanation of some outliers on the bot- 
tom branch may be that they never even proceeded along 
the right and bottom branches to reach their current lo- 
cations. As discussed in § 4.1.1, these may be galaxies 
from the quenched regime that have experienced small 
accretion events, causing them to travel vertically within 
the fueling diagram (i.e., through the bottom branch 
rather than along it) as they develop/deplete small cen- 
tral molecular gas concentrations. In this scenario, such 
galaxies would only fall on the bottom branch as an ac- 
cident of timing. 

4.3. Revisiting the Link Between Mass- Corrected 
Blue- Centeredness and Galaxy Interactions vs. 
Bars 

Throughout this paper, we have made the assumption 
that mass-corrected blue-centeredness is linked to galaxy 
interactions. This assumption is motivated by Kannap- 
pan et al. (2004) , who link blue-centered galaxies to mor- 
phological peculiarities indicative of galaxy interactions. 
One possible issue with this assumption is that E/SOs 
often lack morphological peculiarities (partially by their 
definition of having smooth light distributions), but may 
still have blue centers. Therefore, it is not necessarily ob- 
vious that E/SOs have blue centers for the same reason 



that clearly disturbed galaxies do. 

Our results argue strongly in support of the assump- 
tion that high mass-corrected blue-centeredness implies 
a recent galaxy encounter, even in cases where morpho- 
logical peculiarities are not obvious. In fact, the galaxies 
that have experienced the strongest encounters without 
quenching (gas-rich major mergers of low-mass galaxies) 
are probably the blue-sequence E/SOs on the right and 
bottom branches that show relatively smooth structure 
in their optical images. A notable exception to the link 
between blue-centered color gradients and recent interac- 
tions is the existence of the galaxies in the "dusty" zone 
above the left branch, which appear to be in early, more 
dust-embedded stages of star formation. This implies 
that there is a window shortly after the start of induced 
star formation where mass-corrected blue-centeredness is 
a poor indicator of a recent interaction. 

Bars have been suggested as an alternate mechanism 
for tunneling gas to the centers of galaxies based on di- 
rect observations of gas kinematics (Regan et al. 1995; 
Laine et al. 1999; Regan et al. 1999). The relative im- 
portance of bars versus interactions is not well known due 
to a lack of large, homogeneous samples capable of ade- 
quately testing both mechanisms. Recently, Ellison et al. 
(2011) used the abundance of bars and close pairs in a 
large sample from the SDSS to argue that bars induce 
^3.5 times more central star formation than galaxy in- 
teractions, although the authors note that minor interac- 
tions (i.e., pairs with mass ratios larger than 3/1) are not 
considered in their analysis. Whether bars are a mech- 
anism for inward gas transport completely independent 
of galaxy interactions is also unclear. Galaxy interac- 
tions can induce bar formation (Gerin et al. 1990; Miwa 
& Noguchi 1998), but at the same time bars can form 
in stable disks (Ostriker & Peebles 1973; Sellwood 1981), 
and interactions may actually destroy barred galaxies in 
some cases (Berentzen et al. 2003; Casteels et al. 2012). 
Studies of bar fractions likewise give mixed results as to 
whether bars are related to galaxy interactions or not 
(Aguerri et al. 2009; Li et al. 2009; Barway et al. 2011; 
Lee et al. 2012). 

Bars do not appear to play any role in evolution along 
the right and bottom branches of the fueling diagram. In 
fact, they are absent in these branches, except for the re- 
gion very near the junction between the left and bottom 
branches, which we have already argued is a spiral re- 
birth stage long after the inflow event that triggered the 
starburst. This lack of bars in the blue E/SO population 
is consistent with observations of the bar frequency as a 
function of morphology (Nair & Abraham 2010; Barway 
et al. 2011; Lee et al. 2012). One could argue that bars 
did induce the original gas inflow and their absence on 
the right and bottom branches is caused by the build- 
up of central mass concentrations which dissipate bar 
structures (Norman et al. 1996; Shen & Sellwood 2004; 
Athanassoula et al. 2005; Bournaud et al. 2005). How- 
ever, this picture would not explain why all the galaxies 
are E/SOs. More likely, any existing bars were destroyed 
by violent mergers whose remnants populate the right 
and bottom branches of the fueling diagram. 

Bars are found on the left branch, but the bar frac- 
tion shows no smooth trend with either AC"" or with 
H2/HI (although bars are absent at the top-right of the 
left branch), and they arc even quite abundant in the 



18 



Stark et al. 



bottom-left corner of the fueling diagram where galax- 
ies show no sign of gas inflow. This lack of a smooth 
correlation may be explained if bar lifetimes are on the 
order of a few Gyr (Jogee et al. 2004; Bournaud et al. 
2005; Debattista et al. 2006), in which case bars may 
drive gas inflow that boosts central gas concentrations 
and star formation but also remain well after central gas 
concentrations have been depleted and starbursts have 
ceased (Sheth et al. 2005; Wang et al. 2012, but see Ho 
et al. 1997 and Sakamoto et al. 1999 for alternate view- 
points). This possible longevity makes interpreting the 
role bars play on the left branch difficult. However, since 
there are galaxies near the top of the left branch without 
bars, they certainly do not seem to be required to initi- 
ate an inflow event. Bars may increase the strength of 
inflows and induce quicker depletion, but at this point it 
is unclear whether barred and unbarred galaxies behave 
systematically differently on the left branch. 

5. CONCLUSIONS 

Using mass-corrected color gradients and H2/HI ra- 
tios for a sample of galaxies spanning a broad range 
of morphologies, stellar masses, and evolutionary states, 
we have analyzed the relationship between recent central 
star formation enhancements, likely to reflect galaxy in- 
teractions, and global H2/HI ratios and total gas content 
in galaxies. We summarize our main results: 

• The parameter space of global H2/HI and recently 
enhanced central blueness, which we refer to as 
the "fueling diagram," shows a complex relation- 
ship composed of three main branches - the left 
branch, the right branch, and the bottom branch 
- with most of our galaxies falling on the left 
branch. Galaxies in speciflc evolutionary states 
tend to concentrate in certain regions of the dia- 
gram (e.g., dusty, early-stage starbursts), or can 
be represented on one axis of the diagram (e.g., 
quenched systems). Since our sample is not statis- 
tically representative of the galaxy population, we 
cannot estimate the frequency with which galaxies 
fall on each branch. 

• The left branch is composed primarily of star- 
forming spiral galaxies with a wide range of stel- 
lar masses and gas fractions. It follows a positive 
correlation between global H2/HI and recently en- 
hanced central star formation. We interpret this 
correlation as evidence that H2/HI ratios are sys- 
tematically linked to local encounters with other 
galaxies that drive inflows and replenish molecu- 
lar gas reservoirs. Additionally, apparent enhance- 
ment of H2/HI ratio measurements may be caused 
by decreased Xqo in the high surface density gas 
often found in the centers of galaxies that expe- 
rience inflow events. Galaxies on the left branch 
likely evolve in both directions along it before and 
after inflow events. 

• The right and bottom branches are composed al- 
most exclusively of gas-rich blue-sequence E/SO 
galaxies with stellar masses below the bimodality 
scale Mf, and typically also below the gas-richness 
threshold scale M^. Several lines of evidence sug- 
gest these two branches are part of a continuous 



evolutionary sequence of galaxies formed by gas- 
rich mergers of galaxies with roughly equal masses, 
which results in E/SO galaxies that are experienc- 
ing strong central starbursts, depleting their molec- 
ular gas, and then fading back towards the left 
branch. 

• The population of galaxies on the bottom branch 
evolving back towards the left branch shows a gen- 
eral increase in total gas content and displays a 
clear transition from primarily E/SO to primarily 
spiral morphologies. These results strongly sug- 
gest fresh cosmic gas accretion and post-merger 
disk rebuilding in the low mass regime. Our cur- 
rent analysis does not constrain the timescale of 
this regrowth, but this question will be a topic of 
follow-up research. 

• E/SOs above Mt and especially M;, do not obviously 
move along the branches and may instead move 
vertically in the plot, due to minor accretion events 
associated with nuclear fueling of star formation or 
AGN. 

• Barred galaxies are common on the left branch, 
although the presence of a bar shows no clear cor- 
relation with mass-corrected blue-ccnteredness or 
H2/HI. It is unclear whether bars are involved in 
the small inflow events that drive the evolution on 
the left branch, but bars are likely destroyed in the 
mergers that create the right /bottom branches and 
therefore play little role in these galaxies' evolution. 

The fueling diagram presented in this study links the 
amounts of atomic and molecular gas fuel in a galaxy 
with a metric for the events that drive central fueling and 
HI-to-H2 conversion, providing a useful framework for 
understanding how interactions, inflows, and gas accre- 
tion drive the continued growth and evolution of galax- 
ies. The movement of galaxies through the intercon- 
nected sequences of the fueling diagram highlights their 
dynamic evolution. The left branch of the fueling di- 
agram holds the "normal" star-forming galaxies, which 
appear to progress up and down along the left branch 
during and after inflow events, with star formation al- 
ternately concentrated in the center vs. outer disk. The 
right and bottom branches of the fueling diagram hold 
the more dramatically transforming galaxies that have 
likely experienced recent gas-rich mergers and central 
starbursts. These typically low mass spheroids proceed 
along the right and bottom branches until reconnecting 
with the left branch, potentially re-forming disk galax- 
ies along the way. Some galaxies may proceed through 
the branches of the fueling diagram multiple times un- 
til quenching mergers drive them off the plot. Thus the 
interplay of bulge building and disk regrowth is a funda- 
mental process revealed in varying degrees by the distinct 
evolutionary tracks in the fueling diagram. 
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Table 1 

New IRAM 30m CO(l-O) and CO(2-l) Measurements 









CO(l-O) 










CO(2-l) 














Name 


Flux 


cz 


Wso^' 


rms 


Range ^ 


Flux 


cz 


Wso^" 


rms 


Range "" 


R25 


Log Mh2 '^ 


LogMH2,corr "''^ 






Jy km s^^ 


km s^^ 


km s ^ 


mjy 


km s~^ 


Jy km s^^ 


km s^^ 


kms-i 


mJy 


km s~^ 


II 


Mq 


Mq 




UGC439 


38.81±1.42 


5304 


116±3 


32.17 


5208-5395 


41.91±1.32 


5306 


112±3 


32.06 


5221-5383 


34.0 


9.38±0.02 


9.71±0.02 




UGC1154 


19.46±0.68 


7686 


98±20 


11.40 


7596-7937 


21.46±1.16 


7729 


171±9 


20.85 


7591-7889 


24.2 


9.43±0.02 


9.62±0.02 




UGC1155 


8.21±0.63 


3190 


127±9 


14.15 


3087-3275 


11.95±1.63 


3185 


[121±7] 


37.37 


3113-3295 


20.1 


8.26±0.03 


8.39±0.03 




NGC2780 


16.96±0.92 


1992 


203±6 


17.36 


1863-2132 


20.89±1.53 


1973 


79±19 


33.67 


1886-2084 


31.7 


8.27±0.02 


8.51±0.02 




NGC2844 


28.73±1.70 


1523 


223±22 


25.75 


1300-1720 


29.51±1.83 


1492 


[266±14] 


30.92 


1321-1658 


54.9 


8.31±0.03 


8.69±0.03 




NGC3011 


<1.57 






12.76 


1436-1598 


<4.63 






37.60 


1436-1598 


24.8 


<7.06 


<7.13 




IC2520 


65.96±0.88 


1259 


138±4 


15.44 


1122-1437 


100.50±1.28 


1260 


108±2 


24.57 


1136-1399 


26.6 


8.51±0.01 


8.73±0.01 




UGC5378 


17.09±0.88 


4151 


188±13 


15.32 


3993-4310 


26.49±0.95 


4165 


129±8 


18.89 


4022-4264 


38.2 


8.85±0.02 


9.09±0.02 




NGC3213 


8.62±0.69 


1377 


113±6 


16.18 


1291-1464 


8.26±1.92 


1366 


[83±41] 


32.77 


1231-1560 


36.5 


7.74±0.03 


8.03±0.03 




IC2591 


5.34±1.25 


6706 


[246±36] 


18.92 


6479-6899 












39.2 


8.76±0.10 


9.04±0.10 




UGC6003 


8.10±0.75 


5752 


83±8 


15.24 


5630-5860 












19.1 


8.80±0.04 


8.84±0.04 




UGC6104 


<2.96 






18.88 


2805-3067 


<4.91 






31.40 


2805-3067 


50.9 


<7.83 


<8.11 




NGC3633 


81.45±1.49 


2604 


307±4 


22.93 


2406-2813 


147.24±3.30 


2603 


271±11 


49.66 


2406-2831 


41.6 


9.15±0.01 


9.39±0.01 


CO 


IC692 


1.83±0.30 


1192 


62±9 


12.07 


1161-1222 


<2.36 






24.07 


1106-1208 


24.2 


<6.97 


<7.02 


t^ 


UGC6545 


14.38±0.73 


2657 


94±18 


12.77 


2504-2816 


24.05±0.73 


2641 


105±8 


14.13 


2534-2793 


46.3 


8.45±0.02 


8.72±0.02 


^ 

w 


UGC6570 


14.07±0.58 


1621 


96±7 


12.77 


1505-1701 


24.03±0.56 


1621 


96±3 


13.74 


1539-1696 


35.4 


8.10±0.02 


8.19±0.02 


UGC6637 


2.56±0.51 


1875 


[41±45] 


10.35 


1732-1964 


4.00±0.38 


1866 


100±6 


10.05 


1795-1936 


27.0 


7.45±0.09 


7.50±0.09 


H 


UGC6805 


8.09±0.77 


1159 


121±7 


14.07 


1050-1338 


8.34±1.23 


1175 


49±18 


29.70 


1059-1224 


20.3 


7.56±0.04 


7.60±0.04 


H^ 


IC746 


5.84±0.88 


4978 


[31±29] 


17.99 


4861-5091 


4.46±1.32 


4970 


[75±6] 


40.81 


4924-5025 


43.7 


8.56±0.07 


8.81±0.07 


^ 


UGC7020A 


24.82±0.66 


1519 


110±4 


13.60 


1394-1618 


39.23±1.10 


1520 


107±4 


22.13 


1412-1651 


33.8 


8.29±0.01 


8.37±0.01 


UGC7129 


22.69±0.91 


959 


80±4 


22.47 


886-1043 


25.71±0.91 


959 


55±4 


25.20 


896-1021 


41.9 


7.92±0.02 


8.24±0.02 




NGC5173 


8.98±1.22 


2438 


[154±2] 


25.00 


2302-2531 












35.0 


8.23±0.06 


8.34±0.06 




NGC5338 


13.31±0.89 


845 


64±4 


19.92 


752-943 


29.16±1.62 


849 


59±5 


36.53 


725-915 


56.1 


7.19±0.03 


7.38±0.03 




NGC5596 


<2.15 






15.03 


3058-3276 


<3.11 






21.76 


3058-3276 


30.8 


<7.79 


<7.87 




NGC5762 


<2.61 






19.05 


1688-1888 


<2.99 






21.86 


1688-1888 


44.2 


<7.44 


<7.82 




UGC9562 


3.54±0.95 


1256 


[224±17] 


19.61 


1144-1368 


8.62±2.13 


1330 


[95±47] 


39.18 


1105-1388 


27.5 


7.39±0.12 


7.45±0.12 




UGC9562NE <= 


<1.09 






20.61 


1262-1292 


1.48±0.07 


1272 


[8±4] 


47.12 


1262-1292 










UGC9562SW ' 


<0.78 






16.18 


1035-1060 


1.22±0.13 


1049 


[8±2] 


39.59 


1035-1060 










IC1066 


11.18±0.87 


1572 


107±14 


18.32 


1488-1704 


9.43±0.83 


1566 


[85±14] 


21.18 


1505-1652 


46.0 


7.99±0.03 


8.33±0.03 




NGC5874 


10.72±0.59 


3148 


204±9 


10.60 


3011-3306 


9.77±1.21 


3132 


[133±20] 


18.94 


2990-3383 


64.4 


8.48±0.02 


9.00±0.02 




NGC7328 


48.65±1.71 


2793 


202±17 


28.06 


2656-3014 


46.68±1.50 


2803 


158±16 


25.84 


2675-3000 


62.5 


8.98±0.02 


9.36±0.02 




NGC7360 


<2.29 






12.45 


4517-4877 


<3.79 






20.64 


4517-4877 


36.5 


<8.07 


<8.15 




UGC12265N 


13.29±1.10 


5651 


[246±27] 


17.15 


5435-5832 












16.8 


9.00±0.04 


9.03±0.04 




NGC7460 


40.41±0.81 


3183 


156±4 


16.18 


3052-3295 


41.74±1.11 


3177 


137±4 


21.22 


3054-3315 


44.8 


g.oiito.oi 


9.43±0.01 




NGC7537 


35.37±1.93 


2643 


210±20 


32.05 


2500-2847 


40.80±2.18 


2657 


[218±7] 


39.36 


2538-2834 


68.3 


8.77±0.02 


9.14±0.02 





Note. — ^ Brackets denote galaxies with S/N<6 where linewidths become extremely unreliable and should be used with caution. ^ Range of velocities used in 
integration. '^ Masses include factor of 1.4 to account for Helium. "* Beam-corrected H2 mass (see § 2.3.6). '^ Offset from center of UGC9562 by +9", +11. 5" to observe 
polar ring. Upper limit integration on CO(l-O) flux based on integration range for CO(2-l) detection. CO(2-l) measurements done at a resolution of 2.6 kms~^ due to 
extremely small linewidth. ' Offset from center of UGC9562 by —6",— 9.7" to observe polar ring. Upper limit integration on CO(l-O) flux based on integration range for 
CO(2-l) detection. CO(2-l) measurements done at a resolution of 2.6 kms^^ due to extremely small linewidth. 
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Table 2 

New ARO 12m CO(l-O) Measurements 



Name 


Flux 


cz 


W50" 


rms 


Range ^ 


R25 


Log MH2'' 


Log MH2,corr "'"^ 




Jy km s^^ 


kms-i 


im s ^ 


mjy 


km s"-"^ 


// 


Mq 


Mq 


UGC439 


32.37±2.49 


5301 


129±5 


60.80 


5218-5379 


34.0 


9.30±0.03 


9.39±0.03 


UGC1154 


45.52±4.53 


7770 


194±7 


81.05 


7632-7932 


24.2 


9.79±0.04 


9.84±0.04 


UGC1155 


<22.34 






158.21 


3052-3264 


20.1 


<8.69 


<8.72 


NGC2780 


40.38±4.23 


1976 


204±5 


70.20 


1760-2109 


31.7 


8.64±0.05 


8.70±0.05 


UGC4902* 


5.79±2.49 


1633 


67±18 


93.91 


1600-1667 


40.2 


7.83±0.18 


7.85±0.18 


NGC3032* 


114.11±10.79 


1555 


151±5 


258.50 


1474-1642 


60.0 


8.90±0.03 


9.02±0.03 


IC2520 


105.03±6.92 


1219 


147±1S 


124.77 


1072-1368 


26.6 


8.71±0.03 


8.76±0.03 


UGC5744 


41.90±3.64 


3366 


80±18 


88.97 


3270-3431 


22.8 


9.10±0.04 


9.14±0.04 


NGC3419* 


53.66±4.91 


3040 [107±34] 91.01 


2884-3164 


36.9 


9.05±0.04 


9.07±0.04 


UGC6003 


9.40±1.96 


5833 [184±12 44.79 


5741-5925 


19.1 


8.87±0.09 


8.87±0.09 


NGC3633 


147.45±7.49 


2588 


325±9 


119.83 


2401-2776 


41.6 


9.41±0.02 


9.47±0.02 


UGC6545 


33.21±3.95 


2670 


30±26 


104.81 


2564-2701 


46.3 


8.81±0.05 


8.89±0.05 


UGC6570 


17.86±3.02 


1574 


106±7 


77.30 


1504-1651 


35.4 


8.21±0.07 


8.22±0.07 


NGC3773* 


9.14±2.33 


1002 


[82±9] 


79.92 


961-1043 


35.1 


7.05±0.11 


7.07±0.11 


NGC3870* 


7.96±1.68 


747 


69±13 


62.56 


712-782 


31.2 


7.27±0.09 


7.34±0.09 


UGC6805 


13.79±3.16 


1096 


32±33 


85.13 


1071-1204 


20.3 


7.79±0.10 


7.80±0.10 


UGC7020A 


35.23±4.25 


1536 [138±lll 93.26 


1450-1650 


33.8 


8.44±0.05 


8.46±0.05 


UGC7129 


72.99±7.12 


930 


91±24 


154.42 


822-1027 


41.9 


8.43±0.04 


8.52±0.04 


NGC5173 


<20.31 






152.83 


2373-2561 


35.0 


<8.58 


<8.60 


NGC5338 


23.56±3.46 


807 


40±29 


89.07 


751-896 


56.1 


7.44±0.06 


7.48±0.06 


NGC7077 


<14.56 






142.54 


1086-1198 


25.9 


<7.76 


<7.77 


NGC7328 


98.88±8.56 


2763 [148±55l 134.02 


2592-2984 


62.5 


9.29±0.04 


9.41±0.04 


NGC7328E<= 


37.41±6.42 


2930 


[99±8] 


154.74 


2843-3008 








NGC7328W' 


47.30±6.12 


2763 


217±9 


112.62 


2641-2924 








UGC12265N 


21.33±2.95 


5707 


212±8 


53.57 


5570-5862 


16.8 


9.21±0.06 


9.21±0.06 


NGC7460 


61.89±6.70 


3185 [163±11] 136.14 


3080-3313 


44.8 


9.19±0.05 


9.32±0.05 


NGC7537 


51.20±6.51 


2622 [135±51] 118.60 


2543-2833 


68.3 


8.93±0.06 


9.06±0.06 



Note. — ^ Brackets denote galaxies with S/N<6 where linewidths become extremely unreliable and should be 
used with caution. ^ Range of velocities used in integration. '^ Mass includes factor of 1.4 to account for Helium. 
'i Beam-corrected H2 mass (see § 2.3.6). ° Offset from center of NGC7328 by +28.5", +1.2". * Offset from center 
of NGC7328 by -27", -1.9". * Non-NFGS galaxy. 
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Table 3 

Fueling Diagram Catalog Description 



Column Description 



1 Object name 

2 Right Ascension 

3 Declination 

4 Assumed distance to galaxy 

5 Inclination (see §2.3.1) 

6 u — r color 

7 Stellar mass 

8 Uncertainty in stellar mass 

9 Reference for HI mass 

10 Reference for uncorrected H2 mass 

11 H2 upper limit flag 

12 H2 mass after beam correction 

13 Uncertainty in H2 mass after beam correction 

14 Xco derived from B band luminosity 

15 Xco derived from O/H 

16 u — r blue-centeredness 

17 Uncertainty in u — r blue-centeredness 

18 Stellar mass-corrected u — r blue-centeredness 

19 Uncertainty in stellar mass-corrected u — r blue-centeredness 

20 u — g blue-centeredness 

21 Uncertainty in m — g blue-centeredness 

22 Stellar mass-corrected u — g blue-centeredness 

23 Uncertainty in stellar mass-corrected u — g blue-centeredness 

24 g — '>' blue-centeredness 

25 Uncertainty in g — r blue-centeredness 

26 Stellar mass-corrected g — r blue-centeredness 

27 Uncertainty in stellar mass-corrected g — r blue-centeredness 

28 Morphology (see §3.2.1) 

29 Flag to indicate peculiar galaxy (see §2.3.7) 

30 Flag to indicate presence of bar 

31 Reference for bar classification 



